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Chapter 1 

Introduction 



This course is a major revision of a previous one delivered in Chamonix 
Let us start with a few realistic applications. 

1.1 A few elementary examples 

1.1.1 Linearization of the Riccati equation 

The Riccati equation 



is known to be linearizable, but how to retrieve the explicit formula which 
maps it onto a linear equation? One just looks for the first coefficient uq of 
an expansion for u which describes the dependence on the integration constant 
(a simple pole), i. e. a Laurent series in some expansion function ip(x) : u — 
tp (uo + u\tp + . . .). It is given by balancing the two lowest degree terms 
—Uoip'ip^ 2 = a2(uQ(p~ 1 ) 2 , and the linearizing transformation is then simply the 
change of function u —> <p defined by the singular part transformation u = Mq^ -1 



v! = a2{x)u 2 + ai{x)u + ao(x), 0,2 ^ 0, with ' = — 



(1.1) 




(1.2) 



This will be justified at the end of this course, in section 7.3 



1.1.2 A first integral of the Lorenz model 

The Lorenz model of atmospheric circulation 




, — = xy - bz(x - y) 



(1.3) 
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admits for (b,cr,r) = (0, 1/3, arbitrary) the first integral 



--x 4 + —x(y — x) + (z — r + \)x A 



s 4 */ 3 = K, (1.4) 



but can one go further, i. e. can one obtain more first integrals or even explicitly 



integrate? The answer is yes [ 106 1 . By elimination of (y, z), one first builds the 
second order equation for x(t) 

d 2 z_l/dx\ 2 x 3 K 4t/3 

dt 2 ~ x \dtj 4 3x ■ 1 ' 

For K = this equation admits the first integral 

1 /dx N 2 



and the general solution x = (1/(2A)) cosh(t — to). For K ^ this equation for 
x(t) is equivalent, as shown below, to the following equation for X(T) 

TW/ X' 2 X' aX' 2 + B o 5 , n rN , N 

X" = — - — + ^ + 7 X 3 + -, (a, (3, 7 , 5) constant, (1.7) 

which is the third of six irreducible equations discovered between 1900 and 1906 
by Paul Prudent Painleve and his student Bertrand Gambier and, according 
to the theory of Painleve developed in these lectures, the integration is then 
achieved ("parfaite", says Painleve). Two words may not be familiar to the 
reader : "equivalent" and "irreducible" . "Equivalent" means that the transfor- 
mation law from the physical variables (x, t) to the variables (X, T) which satisfy 
(P3) should not result from a good guess, but should be looked for within a pre- 
cise set of transformations (mathematically the homographic transformations 



(p.5|) defined section 3.3) designed so as not to alter the structure of singulari- 



ties (poles, branch points, . . . ). In this case, one finds that the transformation 
x = a(t)X, T = r(t), with a = — e"*/ 3 , r = ce~ t/3 , c 4 = -±K, (1.8) 

maps the equation for x(t) to the equation (P3) for X(T) with the parameter 
values for (P3) a = = 0, 7 = 5=1. 

As to "irreducible" , it means that there exist no transformation, again within 



a precise class (Drach, Umemura, see section 3.5), reducing any of the six (Pn) 
equations either to a linear equation or to a first order equation. Consequently, 
the general solution of (Pn) has no "explicit expression" , it is just defined by 
the equation itself. There is absolutely no difference between defining the "ex- 
ponential function" from the general solution of u' = u and defining the "P3 
function" from the general solution of the equation (P3). 



6 



1.1.3 A reduction of the Boussinesq equation 

The Boussinesq equation of fluid mechanics 

1 



u t t 



u 2 + -u xx 



-0, (1.9) 



in the stationary case where u does not depend on time, reduces to an ordinary 
differential equation (EDO) which admits two first integrals 

u" + 3u 2 + Kix + K 2 = (1.10) 

and, depending on K\, this ODE is equivalent either to the (PI) equation 

u" = 6u 2 +x (1.11) 

or, after one more integration, to an equation introduced by Weierstrass, the 
elliptic equation 

u' 2 = 4u 3 — g 2 u — 93, (52,53) complex constants. (1-12) 
Both equations have a general solution single valued in the whole complex plane. 

1.2 "Solvable" models, "integrable" equations 
and so on 

Two main fields contributed to the recent interest for the Painleve theory. 
The first one is statistical physics. When Ising solved his one-dimensional model 
and found the partition function F = — (1//3) Log(2 cosh(/3J)), [3 = fc^T, the 
result was a posteriori not surprising. But, when Barouch, McCoy and Wu 
[[37) expressed the correlation function of the two-dimensional Ising model with 
a (P3) function, this strongly contributed to revive the interest for these six 
functions, which now appear in any "solvable model" of statistical physics (see 
Di Francesco, this volume). Retrospectively, the cosh function of Ising is a quite 
elementary output of the Painleve theory. 

The second field is that of partial differential equations (PDE), as shown 
by the above Boussinesq example. After the extension of the Fourier transform 
to nonlinear PDEs []57| , called inverse spectral transform (1ST), Ablowitz and 
Segur H noticed a link between those "IST-integrable PDEs" and the theory of 
Painleve, link expressed by Ablowitz, Ramani and Segur |I| as the conjecture : 
"Every ODE obtained by an exact reduction of a nonlinear PDE solvable by 
the 1ST method has the Painleve property" . For more details, see the book by 
Ablowitz and Clarkson Q and Q. 
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1.3 Insufficiency of quadratures, the need for a 
theory 

Let us return to our main subject, the explicit, analytic integration of ODEs. 
An exceptionally clear introduction ad usum Delphini is the Legon d'ouverture 
(Cteuvrcs vol. I p. 199) given by Painleve in 1895 before starting the Legons 
de Stockholm. For centuries, the question of integration has been formulated 
as : find enough first integrals in order to reduce the problem to a sequence of 
quadratures. But even the simple example of the pendulum shows the insuffi- 
ciency of this point of view. Its motion is reducible to a quadrature defined by 
the integral 



giving the time t as a "function" of the position u. However, this elliptic integral 
does not provide the desired result, i. e. the position as a function of time, 
and, worse, nothing ensures the existence of such an expression. This classical 
problem (the inversion of the elliptic integral) could be solved by Abel and Jacobi 
only by going to the complex domain, leading to a unique value u(t,to,u n ) = 
w,o + sn(t — t 0l k). The symbol sn does deserve the name of function (this is 
one of the twelve Jacobi elliptic functions, equivalent to the unique Weierstrass 
function) because, for any complex k, the application t — > sn(t, k) is single 
valued. 

Following an idea of Briot and Bouquet, this led Painleve to remark (Legon 
no. 1) : "Mais l'importance de cette classe d'equations [a solution generale uni- 
forme] apparait mieux encore si on observe que la plupart des transcendantes 
auxiliaires, dont le role est si considerable (fonctions exponcnticllc, clliptiques, 
fuchsienncs, etc), intcgrcnt des equations diffcrcnticllcs algebriques tres simples. 
Les equations differentielles apparaissent done comme la source des transcen- 
dantes uniformes les plus remarquables, susceptibles notamment de contribuer 
a l'integration d'autres equations diffcrcnticllcs dont Pintegrale n'est plus uni- 
forme." 

This is the famous "double interest" of differential equations : one may 
consider them either as the source for defining new functions, or as a class of 
equations to be integrated with the existing functions available. 



1.4 What can "to integrate" mean? The Painle- 
ve property 



Any converging Taylor series defined on some part of the real line, represent- 
ing for instance a solution of an ODE on some interval —R<x<R, defines 
in fact an analytic function inside the disk |x| < R. Therefore, even when their 




(1.13) 
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variables are real, differential equations and their solutions are naturally defined 
in the complex plane. 

To integrate an ODE is to acquire a global knowledge of its general solution, 
not only the local knowledge ensured by the existence theorem of Cauchy. So, 
the most demanding possible definition for the "integrability" of an ODE is the 
single valuedness of its general solution, so as to adapt this solution to any kind 
of initial conditions. Since even linear equations may fail to have this property, 
e.g. 2u' + xu = 0, u = cx~ x / 2 , a more reasonable definition is the following one. 

Definition. The Painleve property (PP) of an ODE is the uniformizability of 
its general solution. 

Following Bureau |51J, we will call stable an equation with the PP. In the 
above example, uniformization is achieved for instance by removing from the 
complex plane any line joining the two branch points and oo. 

1.5 Singularities of ordinary differential equa- 
tions 

"Les fonctions, comme les etres vivants, sont caracterisees par leurs sin- 
gularites" (Paul Montcl). Singularities are responsible for the limitation of the 
domain of validity of Taylor or Laurent expansions, so their study is mandatory. 

There exits a deep difference between the singularities of solutions of differ- 
ential equations according as whether these equations are linear or nonlinear. 
In the linear case, the general solution (GS) has no other singularities than 
those of the coefficients of the equation once solved for the highest derivative. 
These singularities have a location independent of the arbitrary coefficients of 
integration and they are called fixed. 

On the contrary, solutions of nonlinear ODEs may have other singularities, 
then called movable, at locations depending on the arbitrary coefficients. Thus, 
the equations |89(| 



du u 



dx 



+ = 0, u = ce 1/x , (1.14) 



x 



du u 2 n 1 

— + — = 0, u= , (1.15) 

da; a; c + Log x 

2 

= 0, u — sin(c + Logx), (1-16) 



du yl — u 
dx x 



where c is the arbitrary constant of integration, all have a fixed singularity in 
their general solution at x — (isolated essential singular point for the first one, 
logarithmic branch point for the two others). In addition, among the last two 
ones, which are nonlinear, the second one has movable simple poles, and the 
third one has no movable singularity. All three have the PP. 
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The possible singularities of differential equations have been classified by 
Mittag-Leffler : in addition to the familiar ones (poles, branch points, essential 
singular points), there can exist essential singular lines, analytic or not, or per- 
fect sets of singular points, as illustrated by the Fuchsian and Kleinean functions 
of Poincare. One example is Chazy's equation of class III 

v!" - 2m" + 3m' 2 = 0, (1.17) 

whose general solution is only defined inside or outside a circle characterized by 
the three initial conditions (two for the center, one for the radius); this solution 
is holomorphic in its domain of definition and cannot be analytically continued 
beyond it. This equation therefore has the PP, and the only singularity is a 
movable analytic essential singular line which is a natural boundary. 

1.6 Outline and basic references 

For the outline, we refer the reader to the detailed table of contents; the 
choice made is to develop the construction of necessary conditions, at the ex- 
pense of the explicit integration methods, only briefly introduced in chapter 

The basic texts are due to Painleve and his students and we abbreviate their 
references as Legons ( |89|] L egons de Stockholm, 1895, a high level course on non- 
linear ODEs), BSMF (]90[ first memoir, 1900, on first and second order ODEs), 
Acta ([|2| second memoir, 1902, on second and higher order ODEs), CRAS (Q, 
1906, an addendum after Gambier discovered the functions (P4), (P5), (P6)), 
Gambier (@ these, 1909, on second order ODEs), Chazy (|§ these, 1910, 
on third and higher order ODEs), Gamier (]58| these, 1911, on higher order 
ODEs). Most Painleve works are reprinted in Oeuvres (|9^] three volumes 1973, 
74, 76, again available from CNRS-Editions, e-mail editions@edition.cnrs.fr). 
For a global overview of these results, see the book of Hille [|33"| . preferably to 
the one of Ince |6(|. For a detailed exposition (indeed, in the classical period, 
it was kind of fashionable to avoid details) and additional results, see the three 
memoirs of Bureau M. I [pi M. II |]|, M. Ill jl(|. Peter Clarkson maintains 
an extensive bibliography [|26[ covering both the classical and the recent period, 
reproduced in fl. 
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Chapter 2 

The meromorphy 
assumption 



2.1 Specificity of the elliptic function 

A very deep result of L. Fuchs, Poincare cf. Qeuvres de Painleve III 

p. 189) and Painleve (Legon no. 7 p. 107) is that the class of first order ODEs 
F(u',u,x) = 0, with F polynomial in u' and u, analytic in x, defines one and 
only one function, from the general solution of ( 1.1 2| ) . This function is not 



historically new since this is precisely the elliptic function p introduced earlier 
by Weierstrass, i. e. the particular solution of 

p' 2 =4p 3 - 52 p-ff3 = 4(p-ei)(p-e 2 )(p-e 3 ), (fl-2, 33, e a ) G C, (2.1) 
which admits a pole at the origin 

p(x, g 2 ,g 3 ) = x- 2 + 9 ^x 2 + gx 4 + 0(x e ). (2.2) 

The novelty of p is elsewhere : this is the transcendental dependence of 
the general solution p(x — xq, <?2, 33) on the arbitrary constant xq, which makes 
impossible the reducibility of the elliptic equation to a linear equation. Among 
the many nice properties of elliptic functions (see e.g. Q), the most interesting 
to us is their structure of singularities. These are doubly periodic meromorphic 
functions (which is their usual definition), and there exists an entire function 
tr, i. e. without any singularity at a finite distance, whose — p is the second 
logarithmic derivative 

P = C = Log a, C" 2 + 4C' 3 - 92C + ff 3 = 0. (2.3) 

ax ax 
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Therefore the only singularities of the general solution p(x — xq, 93) of (1.12) 
come from the zeroes of a and are a lattice of movable double poles located at 
xq + 2mu> + 2nui', with m and n integers, u>, u> the two half-periods. 



2.2 The meromorphy assumption 

The Laurent expansion ( [2.2] ) certainly motivated two students of Weierstrass, 
Paul Hoyer and Sophie Kowalevski, to investigate further the possibility for the 
GS of an ODE to be represented by a Laurent series with a finite principal part, 
so as to exclude essential singularities. This meromorphy assumption, briefly 
said, consists in checking the existence of the Laurent series and its ability 
to represent the GS, i. e. to contain enough arbitrary parameters. But, since a 
Laurent series is only defined inside its annulus of convergence, this study is only 
local and it cannot dispense of a further study in order to explicitly integrate, 
using completely different means. 

The first attempt is due to Hoyer in 1879 164] with the system 









-( 







ax 



Cl 




(2.4) 



under the restriction that neither the determinant nor any of its first or sec- 
ond order diagonal minors vanishes; he even generalized the assumption to the 
Puiscux series 



E 

3=0 



A tJ {(t-t )-}- n +i, i= 1,2,3, 



with n and r positive integers and (j4io, ^20j ^30) 7^ (0,0,0), but in fact the 
numerous cases of integrability by elliptic functions which he discovered were 
found by a direct Ansatz, and not as necessary conditions for the Laurent series 
to exist. Continued by Kowalevski with a quite similar system except that it 



is six-dimensional, which will be seen section 2.5, the method will only get its 



final shape with Gambier in 1910 (see pages 9 and 49 of his thesis). 

2.3 A flavor of the meromorphy test 

We must warn the reader that this section is not the algorithm to apply, but 



just a flavor of it; the final algorithm will only be given section 3.£. 

Let us start with a single equation; the case of a system is not different, 
apart from technical complications. 



Consider the equation 



d 2 u 



E{x,u) = -^ + 6u 2 +g{x) = 0, (2.5) 



12 



with g analytic. 

Assume that u has a polar behavior at some location xo distinct from any 
of the possible singularities of the coefficients of the equation, here g{x); such a 
pole is therefore movable. One has to check the existence of all possible Laurent 
series with a finite principal part 

+00 

u = ^2 u jX 3+p , X = x ~ x o, u ^ 0, (2.6) 

in which — p is the order of the pole, which must be an integer, and the coeffi- 
cients Uj are independent of x. 



After insertion of this series in (2.5), which is polynomial in u and its deriva- 
tives, and replacement of g{x) by its Taylor series in the neighborhood of Xq, 
the left-hand side, as a sum of Laurent series, is itself a Laurent series with a 
finite principal part 

E = [p(p-l)u oX P- 2 + (p+l)pu 1 x p - 1 + ...}+6[u 2 x 2p + 2u u lX 2p+1 + ...} 
+{g(x )+g'(x )x +•••], (2-7) 

which we denote more generally 

+00 

E = Y,EiX j+ \ (2.8) 

q being the smallest integer of the list (p — 2,2p, 0). The method consists of 
expressing the conditions for this series to identically vanish : Vj € N : Ej = 0. 



First step. Determine all possible families of movable singularities (uo,p). 
This is expressed with three conditions : 

• 1) (condition u ^ 0) equality of at least two elements of the list (p — 
2,2p, 0) (q denotes their common value), the involved terms of E being 
called dominant and denoted E, 

• 2) (dominance condition) inferiority of q to the other elements of the list, 

• 3) (vanishing Laurent series condition) vanishing of the coefficient Eq of 
the lowest power \ q , which involves only the dominant terms 

E = 15m X ~ q E{x, u x p ) = 0, u + 0, (2.9) 



i. e. respectively : one linear equation for p by pair of terms considered, several 
linear inequations for p, one algebraic equation for (uo,p). 
A necessary condition to prevent multivaluedness is then 
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• CO. All possible values for p are integer. 

If there exists no family which is truly singular (p negative), the method 
stops without concluding. 

Here, the various possibilities for these linear equations and inequations are 

q = p-2 = 2p and q < 0, (2.10) 
q = 2p = Q and<j<p-2, (2.11) 
q=p-2 = {) and q < 2p. (2.12) 

Their geometric representation is known as the Puiseux diagram or Newton's 
polygon (ref. Q sec. 3.3, |36| sec. 12.61). The two solutions (p,q) = (—2,-4), 
(2, 0) satisfy the condition CO and the second one must be rejected, as being 
nonsingular. So, the dominant part is here E = —u" + 6it 2 . The algebraic 
equation Eq = 

E = -6u + 6u 2 = 0, u ^ 0, (2.13) 

has only one root uo = 1. 

For j = 1,2,.. ., each successive equation Ej — has then the form 



Vj > 1 : Ej = P(u ,j) Uj + Q^m | I < j}) = 0, (2.14) 



here 



P(uo,j) = ~U - 2)(j - 3) + 12uo - -(j + l)(j - 6), (2.15) 
Qi =0, Q 2 = 6u?, Q 3 = 12ui«2, (2.16) 

Vi>4: Q^^M + e^ufett^fc. (2.17) 
^ fe=i 

So the sequence Ej = 0,j > 1, is just one linear equation with different right- 
hand sides, and it can be solved recursively for uj. Whenever the positive integer 
j is a zero of P, two subcases occur : either Qj does not vanish and the Laurent 
series does not exist, or Qj vanishes and Uj is arbitrary. Since Xq is already 
arbitrary, in order to represent the GS, one wants N — 1 additional arbitraries 
to enter the expansion, where N is the order of the ODE. Let us admit for a 
moment that the value j — —1 is always a zero of P, a result whose general 
proof (given section 5.5) needs prerequisite notions of perturbation theory; this 
value j = — 1 will be seen to represent the arbitrary location of xq. Hence the 
following steps. 

Second step. For each family, determine the polynomial P (do not compute 
Qj yet) and require the necessary conditions : 

• CI. The polynomial P has degree N. 

• C2. N — 1 zeroes of P are positive integers. 
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• C3. The N zeroes of P are simple (i. e. of multiplicity one). 



If either CI, C2 or C3 is violated, the method stops and one concludes to 
a failure because the general solution cannot be meromorphic. 

The zeroes of P are called indices and P — itself is the indicial equation. 



Indeed, anticipating on the exposition of the general theory sections 5.6 and 
|5.7| , they are the Fuchs indices i near x = of a linear equation introduced by 
Darboux fl5| under the name "equation auxiliaire" , so the indicial equation is 
computed as follows |52|] . Take the derivative of E(x,u) with respect to u 

. , E(x. u + Xv) — E(x, u) 
Vv : E'(x,u)v ee lim , ~ J ~ J -, (2.18) 

A — >0 A 

here 

E'{x,u) ee -dl + 12u; (2.19) 

evaluate this linear operator at point u = UoX p , which defines the "auxiliary 
equation" (i. e. the linearized equation at the leading term) 

Vw : E'(x,u oX p )v = 0, (2.20) 

here 

Vt- :E'(x,u oX - 2 )v = (-dl + l2u oX - 2 )v = 0; (2.21) 
establish the indicial equation of this linear ODE near its Fuchsian singularity 

x = o 

Pii) = lim x~ l - q E'(x, u oX p )x l+p = 0, (2.22) 

here 

P(i) = limx- l+ H-d 2 x + l2uo X - 2 )x 1 - 2 
= -(i-2)(i-3) + 12uo 

= — (i + l)(i — 6). (2.23) 

The shift i — > i + p in the above equation is just a convention aimed at not 
producing an unfortunate difference between the Fuchs index i and the index j 
of the recursion relation Ej =0. 

Now, one just has to check the existence of the Laurent series. 

Third step. For every positive integer zero i of P (a Fuchs index), require 
the condition 

• C4. 

Vt € TV, P(i) = :Qi = 0. (2.24) 
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This is done by successively solving the recursion relation up to the greatest 
positive integer Fuchs index. As soon as a C4 condition is violated, one stops 
and concludes to a failure : the ODE has not the PP. After the greatest positive 
integer Fuchs index has been checked, the method is finished. 

Here, one finds 

u = 1, ux = u 2 = u 3 = 0, it 4 = -7^, u 5 = -^r, (2.25) 

10 

and the condition C4 at index i = 6 is 

Q 6 ee g»/2 = 0. (2.26) 



i.e., since xo is arbitrary, g" = 0. The ODE (|J) is restricted to be fll.lOD 

which 

has been seen to have a meromorphic GS, so in this case the generated necessary 
conditions are sufficient. 
Remarks. 

1. We have retained the classical vocabulary ("famille" is used by Gambier, 
p. 38 of his thesis Q, "indices" is used by Gambier, Chaz y Q and Bu- 
reau [[l3)), rather than the one more recently introduced j^, |4| ( "branch" , 
"resonances"). Indeed, "branch" has another meaning in classical analysis, 
where it denotes a determination of a multivalued application, which may 
create some confusion. As to "resonance", its identification with a basic 
notion of a linear theory, the Fuchs indices, makes useless the introduction 
of such a term. 

2. Conditions Qi = at Fuchs indices i are often referred to as "no-logarithm 
conditions" because, if some of them are not satisfied, there exists a gen- 
eralization of the Laurent series, called ip— series (||| chap. 7), which is 
a double expansion in \ an( i Log%. This series contains no logarithms 
(i. e. reduces to the Laurent series) iff all Qi vanish. 

3. One must prove that the radius R of the punctured disk \x — xo\ < R in 
which the series converges is nonzero. 

4. As indicated by Gambier |5(| p. 50, there is no need to expand the coeffi- 
cients g(x) of the equation around Xo- This is achieved fl30| by taking for 
the expansion variable not x — Xq, but a mute variable x with the only 



property \x — 1- Coefficients Uj in eq. (2.25) become dependent on x 
instead of xq : 

1 n 9{x) g'{x) 

Uq = 1, Ul =U 2 = 1*3 = 0, Ui = — — , u 5 = — - , 

10 15 
Q 6 ee ^ = 0. (2.27) 
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Let us insist again on the danger of using the present test as it is. An 
example of Chazy makes evident the necessity for a more reliable test : the 
equation with a single valued general solution (|22| p. 360) 

{v!" ~ 2u'u") 2 + Au" 2 {u" - u' 2 -l) = Q, u = e Cl *+ c Vci + ^f^x + c 3 . 

possesses a logarithmic family u ~ — Log(a; — Xq). 
Exercise 2.1 Handle the equation 

2uu" - 3u' 2 = 0, u = cx/{x - c 2 f. (2.28) 

Solution. 

2p — 2 = 2p — 2, Eq = 2uqp(p — 1) — 3u^p 2 — 0. Hence p = —2, uq arbitrary.D 

(2.29) 



2.4 Extension to a system 

If the differential equation is defined by a system 

E(a;,u) = 0, (2.30) 

(boldface characters represent multicomponent quantities) , the scalar equations 
of section 2.3 become systems : a linear system for the components of p, an 
algebraic system for the components of Uo, a linear system with a rhs for uj, a 
determinant for the indicial equation, etc. Take the example of the Euler system 
(diagonal Hoyer system) 

da;i dx2 dx^ 

Ei = -r--ax 2 x 3 = 0, E 2 = ——-px 3 xi = 0, E 3 = —— -7x1X2 = 0. (2.31) 
at at at 

First step. The necessary condition on (p, uo) is 

• CO. All components of p are integer, all components of uo are nonzero. 

Of course, if a component of Uo is zero, one must increase by one the as- 
sociated component of p until the new component of Uo becomes nonzero. It 
there exists no truly singular family (at least one component of p negative), the 
method stops without concluding. 

Here, the unique solution p of the linear system qi = p\ — 1 = p 2 + pj, and 
cyclically is thus p\ = p 2 = p 3 = — 1, q\ = q 2 = q 3 = — 2. The algebraic system 
for Uo 

E = lim x" q E(x, u oX p ) = 0, u ^ 0, (2.32) 
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is written 

Ei,o = — xi,o — ax2flx?,fl = and cyclically, (2.33) 
and it defines four families 

x i,o = x l,o = — - ^3.0 = 4s i a;i,oa;2,oa;3,o = (2.34) 

P7 ' -fa ap apf 

which we gather under the unique algebraic writing 0:1,0 = a, X2.0 — b, 2:3,0 = c. 
Second step. The linear system 

Vj > 1 : Ej = P(j) Uj + Q,({u ; I i < i}) = (2.35) 

generates the indicial equation 

dct P(») = 0, P(i) = lim x-'- q E'(a;, u aX p )x l+p , (2.36) 



i — 1 -ac — ab \ 

-/3c i-1 -/3a , detP(i) = (i + l)(i-2) 2 = 0. (2.37) 
—76 —7a i — 1 J 

Classical results from linear algebra on the resolution of the matrix equation 
AX = B give the necessary conditions 

• CI. The polynomial detP has degree N. 

• C2. N — 1 zeroes of det P are positive integers. 

• C3. Every positive zero i of det P has a multiplicity equal to the dimension 
of the kernel of det P(i). 

Here, each of the four families has the same indices (—1,2,2) and, for the 
double index i = 2, the three rows of matrix P(2) are proportional, so its kernel 
has dimension two. 

Third step. For every positive integer zero i of P (a Fuchs index), require 
the condition 

• C4. 

Vi € A/", dctP(i) = : the vector is orthogonal to the kernel 

of the adjoint of operator P(i). (2.38) 

Here, the condition C4 is satisfied at index two, and the Laurent series are 
finally 

Xl =a X - 1 +a 2X + 0( X 2 ), X = t-to, (2.39) 
x 2 = b X - 1 +b 2X + 0(x 2 ), (2.40) 
x 3 - c X - x + c 2X + 0( X 2 ), a 2 + b 2 + c 2 = 0, (2.41) 

with (t , b 2 ,c 2 ) arbitrary. 
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2.5 Motion of a rigid body around a fixed point 



It is ruled by the system 

B ^T + ( A ~ CVawi + {xik 3 - x 3 h) = 0, ^ - ^ifc 3 + w 3 fci =(0,42) 

dt 

C~7T + (-B - ^4)^1^2 + (xzki - xik 2 ) = 0, - u> 2 ki + u)\k 2 = 0, 
dt at 

depending on six parameters : the components (A,B,C), positive, of the diag- 
onal inertia momentum / and the components (x±, x%, x 3 ), real, of the vector 
linking the fixed point O to the center of mass G. Because it admits the 
three first integrals 

K x = (iTt).Ti - 20G.~t = Auj\ + Bwl + CloI - 2(x 1 k 1 + x 2 k 2 + x 3 k 3 ), 
K 2 = (Iit).'t = Awiki + Buj 2 k 2 + Cuj 3 k 3 , 
K 3 = It Hi = k\ + k\ + k 2 , 

and a last Jacobi multiplier equal to 1, 

' dujj \ v-^ „ / dk 



v dt J ^ \ dt 

a sufficient condition of reducibility to quadratures (i. e. to separation of vari- 
ables, which implies neither meromorphy nor single valuedness) is the existence 
of a single additional first integral independent of time. 
Before Kowalevski, the only such known cases were 

• the isotropy case, with 

A = B = C : K i = O&.Tt = XxU)x + x 2 u 2 + x 3 uj 3 , (2.44) 

• the case of Euler (1750) and Poinsot(1851) G at the fixed point O (x\ = 
x 2 — x 3 = 0) with 

G = : K± = \lTt\ 2 =A 2 uf + B 2 u 2 1 + C 2 ul (2.45) 

• the case of Lagrange (1788) and Poisson (1813) 

A = B, xi = x 2 = : K 4 = u 3 , (2.46) 

and for these three cases the general solution is elliptic, hence meromorphic fl59|. 
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Let us denote a family as (t — to is abbreviated as t) 

+00 +00 
= } j u>i > jt ni+J , h = ^ki,jt mi+3 , wi :0 ^2 : oW3 : ofci,ofe,ofc3,o 7^ 0, 

3=0 j=0 

with Z = 1, 2, 3 and {wij, hj) complex. There exist numerous families, some of 
them with n; — ti/^to/ — not integer. To shorten, let us restrict to the case 
where all the differences ni ~ nk,mi ~ rrik are integer, and redefine a family as 

+00 

ui = t n ^2ui tj t\ (wi,o,w 2 , 0,^3,0) ^ (0,0,0), 
3=0 

fc ; =t m ^fc ; ,^, (^,0^2,0,^3,0)^(0,0,0). (2.47) 
i=o 

One such family is defined |72| by the exponents n/ = —l,mi — —2, and the 
sextuplets (w;,o, fc;,o) solutions of the algebraic system 

AuJlq + (B - C)u) 2fl UJ 3fi + ^2^3,0 - 2^2,0 = 0, 2fci, + W 3i0 fc 2 ,0 ~ ^2,0^3,0 = 0, 

Bu 2 fl + (C - A)u 3fl u lfi + x 3 k lfi - xik 3 . = 0, 2fc 2 , + ui,ok 3 ,o ~ ^3,0^1,0 = 0, 
Cu 3fi + (A - B)u lfl u 2 fi + xik 2 . - x 2 ki. = 0, 2k 3 . + u 2 ,ok lfi - ^1,0^2,0 = 0, 



and the linear system for j > 1 is 



(J 


-l)A 


[C- 


- B)u 3fi 


(C- 


- B)lu 2 . 





%3 


-x 2 \ 




{A- 


- C)w 3fi 


(j 


-l)B 


(A- 


- C)u)i t o 


-x 3 





Xi 


W 2j 


(B- 


- A)ui 2fi 


(B- 


- A)u lfi 


(J 


-1)C 


x 2 


-xi 





w 3j 









k 3} o 




-k 2 , 


J -2 


—^3,0 


^2,0 






-&3,o 









h,o 


^3,0 


J -2 


-^1,0 


k'2,j 

V fc 3 j / 




k2,0 




-ki,o 







^^2,0 


Wl,0 


3-2j 




= 0. 
















(2.48) 



The determinant det P must have five positive zeroes. 

In the generic case (^4, B, C) all different and G ^ O, there exists a unique 
solution to the algebraic system, depending on one arbitrary parameter and the 
root of an eighth degree equation |73|], but the determinant 

detP = ABCU + - 2)(j - 4)(j 2 - j - M ), (2.49) 

where /i is an algebraic expression of (A, B,C,xi,x 2 ,x 3 ), has five positive integer 
zeroes iff /i = 0, which corresponds to inadmissible values for the six parameters 
{A,B,C must be real positive, x\,x 2 ,x 3 real). 

A thorough discussion of the nongeneric cases of this family rt; = — l,mj 
2 led Kowalevski to retrieve the three known cases, as expected, and finally 
to find the subcase 

A = B, (x u x 2 ) ± (0, 0), ulo + c4 - 0, (2.50) 
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for which the unique solution is 



iC . „. . 2 

2{Xi + 1X2) A 

2C 

fei,o = — : — , ^2.0 = «fel,0, ^3.0 = 0, (2-51) 

X\ + 1X2 

det P = ABC(J + l)(j - 2)(j - 3)(j - 4)(j + 1 - 2C/A)(j - 2 + 2C/A), 

in which A is defined by the relation 

2C-A- 4A.t 3 = 0, A ^ 0. (2.52) 

There exist five positive integer indices iff A = 2C, x 3 = 0, and the first integral 

A = B = 2C, x 3 = : K i= \C(u) X + iuj 2 ) 2 + {x x + ix 2 ){k 1 + ik 2 )\ 2 (2.53) 

terminates the proof of reducibility to quadratures. Sophie Kowalevski then 
managed to explicitly integrate with hyperelliptic integrals and to prove the 
meromorphy of the general solution, a feat which won her an instantaneous 
fame. 

Remark. Neither Hoyer nor Kowalevski enforced conditions C3 and C4. 
This was done for the first time by Appelrot |J, who found another family 
ni = —1, mi = with the indices (—1, —1, 0, 1, 2, 2) and, despite the absence of 
five positive integer indices, computed the next terms and found at the double 
index 2 the no-log condition 



x x ^JA{C -B)+ x 2 \J B(A - C) + x 3 y/C(B - A) = 0, (2.54) 

whose real and imaginary parts yield, with the convention A > B > C, 

x 2 = 0, x x ^jA(B -C)+ x 3y /C(A -B)=0. (2.55) 

Nekrasov [^6| and Lyapunov then proved the multivaluedness of this case by 
exhibiting yet another family with complex exponents. 



2.6 Insufficiency of the meromorphy 

Here is the opinion of Painleve (|^2) pp. 10, 83, Oeuvres III pp. 196, 269) : 
"M c Kowalevski se propose de trouver tous les cas ou le mouvement du solide est 
dcfini par des fonctions meromorphes de t qui possedent effectivement des poles. 
Son procede laisse echapper les cas ou ces fonctions seraient uniformes sans avoir 
de poles, soit qu'elles fussent holomorphes, soit que toutes leurs singularites 
fussent transcendantes. 

De plus, apres avoir forme les conditions pour qu'il existe des poles mo- 
biles, M c Kowalevski remarque que ces conditions entrainent l'integrabilite des 
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equations du mouvement, ce qui lui permet de mener la question jusqu'au bout. 
Mais cette remarque laisse echapper un cas oil il existe des poles et qui n'est pas 
un cas d'integration. Toutefois les geometres Russes ont montre, par la suite, 
que, dans ce cas, les equations du mouvement n'ont pas leur integrate uniforme. 

Les resultats de M e Kowalevski subsistent done en fait. Mais, si interessante 
que soit la voie suivie par M e Kowalevski, il etait desirable de reprendre la 
question d'unc facon plus rationncllc. C'cst ce que permettent les precedes 
que j'ai employes pour les equations du second ordre : ils fournisscnt de la 
maniere la plus naturelle ct la plus simple les conditions necessaires pour que 
ce mouvement soit represents par des fonctions uniformes de t, sans qu'il soit 
besoin de faire aucune hypothese sur ces fonctions. Les conditions auxqucllcs 
on parvient ainsi nc different pas d'aillcurs de ccllcs de M c Kowalevski. Pour ce 
probleme particulier, on n'arrive done pas a des cas nouveaux." 

Painleve thus insists that poles are not privileged : they are just one kind 
of singularity among many possible others. 

"Une discussion qui ecarterait d'avance certaines singularites comme in- 
vraisemblables serait inexistante." (Painleve, |92j p. 6, Qsuvres III p. 192). 

"In the statement of the problem, poles are not mentioned; if in the final 
result the particular integrals prove to be meromorphic, it is a result of the 
research. Likewise, no mention is made of one or another type of critical or 
singular point." (Bureau jl9| p. 105). 

Thus, definitely, the meromorphy assumption has to be waived as a global 
property, although it may be, and indeed is, quite useful at the local level. The 
only relevant property based on singularities is the Painleve property as defined 
in section 1.4, and the goal is to build a rigorous theory without any a priori 



on the movable singularities. Despite the pessimistic opinion of Picard who 
thought the task impossible, Painleve built that theory by a clever application 
of the theorem of perturbations of Poincare and Lyapunov. 



2.7 A few examples to be settled 

The theory to come and the resulting "Painleve test" should be able to han- 
dle the following differential systems for which the meromorphy test of section 
[2.3| is inconclusive or even erroneous. We give the location where the solution 
can be found. 

1. ("Soit qu'elles fussent holomorphes^") Extend the test to handle the equa- 
tion 2uu" — u' 2 — 0, with general solution u — (c 1 x + c 2 ) 2 . Solution section 

2. ("Soit que toutcs leurs singularites fussent transcendantes"). Extend the 
test to handle the equations uu" — u' 2 = and 2u 2 u'u"' — "iu 2 u" 2 + u' A = 0, 
whose general solution is, respectively, u = e ClX+C2 and u = cie 1 ^ C2X+C3 \ 
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One may notice that the second equation is the Schwarzian derivative of 
Log it. 



3. The "uncoupled" system with a meromorphic general solution 



du 9 „ dv 
— +u 2 = 0, — 
dx dx 



v 2 = 0. 



(2.56) 



admits two families (modulo the exchange of u and v) 



(Fl) 
(F2) 



u ~ x 
u ~ x 



indices (-1,-1) (2.57) 
v oX°> v o arbitrary, indices (— 1, 0), (2.58) 



of which the first one fails the condition C2. Solution section p. 4 . 

("Un cas ou il existe des poles et qui n'est pas un cas d'integration"). The 
Bianchi IX cosmological model 



(Log A)" = A 2 -(B- C) 2 and cyclically, ' = d/dr, 



(2.59) 



admits for B = C a particular four-parameter meromorphic solution [ 108 1 
k\ n ^ A;| sinh fci (t — Ti ) 



A 



sinhfci(r — n) 



, B = C 



k\ sinh k%{r — T2) 



(2.60) 



Prove the absence of the PP by studying the family p = ( 0, —2, —2) which 
has only four Fuchs indices (— 1, 0, 1, 2). Solution section 5.8.3 . 



5. The Bianchi IX model ( 2.59 ) admits a family p = (—1,-1,-1) with the 
indices (—1, —1, —1, 2, 2, 2). P rove the absence of the PP by studying this 
family. Solution section 5.7.6. 



6. The Chazy's equation of class III (1.17) admits a Laurent series which 
terminates u = Uq/(x — xq) 2 — 6/(x — xq), with (xo, wo) arbitrary. From 
the study of this family, decide about the meromorphy of the general 
solution. Solution section 5.8.1. 

7. In a problem in geometry of surfaces, Darboux Mm encountered the system 



dxi/dt = X2X3 — x\(xi + 2; 3) and cyclically, 



(2.61) 



explicitly excluded by Hoyer, cf. section 2/2, and found the two-parameter 
meromorphic solution 

xx =c/(t-t ) 2 + l/(t-t ),x 2 =x 3 - l/(*-*o), (*o,c) arbitrary. (2.62) 

For c = 0, the Fuchs indices are (—1, —1, —1). Extend the test to build 



no-log conditions at this triple —1 index. Solution section 5.7 
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The equations 



-2uu" + 3m' 2 + d 3 u 3 = 0, d 3 ^ 0, 
u'" + uu" - 2u' 2 = 0, 
u"" + 2uu" - 3u' 2 = 0, 



(2.63) 
(2.64) 
(2.65) 



have no dominant behavio ur. Prove the absence of the PP for each of 
them IE4||. Solution section 5.G. 
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Chapter 3 

The true problems 



In this chapter, we state the true problems and manage logically to the only 
correct definition for the Painleve property (PP) : "absence of movable critical 



points in the general solution", equivalent to that already given in section 1.4. 
This includes 

- the two classifications of singularities of differential equations (fixed or mov- 

able, critical or noncritical) , 

- the two differences between linear and nonlinear (movable singularity, singu- 

lar solution), 

- the statement of the ambitious program proposed by Painleve, a first, quick 

look at the method of resolution (the "double method" and the "double 
interest") and the results for first order (equation of Riccati, function of 
Weierstrass) and second order (classification of Gambier, the six Painleve 
functions). 

All the ODEs considered are defined on C or on the Riemann sphere (i.e. the 
complex plane compactified by addition of the unique point at infinity). 
Firstly, a more precise definition of the term "to integrate" is required. 

Definition. To integrate an ODE, in the "modern sense" advocated by 
Painleve, is to find for the general solution a finite expression, possibly multi- 
valued, in a finite number of functions, valid in the whole domain of definition. 

The important terms in this definition are "finite" and "function" . 

Example 1 (nonintegrated ODE) : the Taylor series u = uq Ylj^o i~ ( x ~ 
fo) u o]"^ f° r the Cauchy solution does represent the general solution of the ODE 
v! + u 2 = but this representation is local and the integration cannot be con- 
sidered as achieved until one has : found the radius of convergence, performed 
the summation, analytically continued the sum everywhere this is possible, and 
identified the analytic continuation with the meromorphic function (x — Xi) . 
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Example 2 (integrated ODE) : the ODE 2uu' —1 = has for general solution 
u = (x — xo) 1 / 2 , a multivalued finite expression built from the "multivalued 
function" (see below) z — > z 1 / 2 . 

Representations by an integral, a series or an infinite product are acceptable 
iff they amount to a global, as opposed to local, knowledge of the solution. 

A prerequisite to the integration in the sense of the above definition is there- 
fore to extend the set of available functions, to serve as a reservoir from which to 
build finite expressions. At this stage, one must go back to the term "function". 

Definition A function is an application of a set of objects into a set of 
images which applies a given object onto one and only one image. 

In other words, a function is characterized by its single valuedness, and terms 
such as "multivalued function" should be carefully avoided. In our context, a 
function is a single valued application of the Riemann sphere onto itself. 

Definition (Painleve [ pp| p. 206). A critical point of an application of the 
Riemann sphere onto itself is any singular point, isolated or not, around which 
at least two determinations are permuted. Common synonyms are : for critical 
point, branch point, point of ramification; for determination, branch. Such a 
point is an obstacle for an application to be a function. 

Examples : the applications x — > \fx — a and x — > Log(a; — a) both have 
exactly two critical points, a and oo. Around each of them are permuted re- 
spectively two determinations and a countable infinity of determinations. 

Remark. An essential singular point is not necessarily a critical point, since 
essential singularities, isolated or not, can be critical or not. Examples of critical 
essential singularities are : x = for tan(Loga;) (nonisolated) or sin(C + Logx) 
(isolated and transcendental, Legons de Stockholm pp. 5-6, |9^ |, |6(| §14.1 
p. 317). Examples of noncritical essential singularities are : x — oo for e x 
or equivalently x = xq for e 1 /^-^ (isolated), x = oo for tanx (nonisolated). 
Although, according to a classical theorem of Picard, an analytic function can 
take any value but at most two (oo and another one) in the neighborhood of an 
isolated essential singularity, a noncritical essential singularity is not an obstacle 
to single valuedness. 

3.1 First classification of singularities, uniformi- 
zation 

Definition. The first classification of singularities is the distinction critical 
or noncritical between singular points of applications. Note that it does not 
involve differential equations. 

Consider a multivalued application of the Riemann sphere onto itself. There 
exist two classical methods, called uniformizations, to define from it a single 
valued application, i.e. a function. 
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The first one is to restrict the object space by subtracting some lines, called 
cuts, so as to forbid local turns around critical points; for the above two exam- 
ples, one removes any line joining the two points a and oo. 

The second method is to extend the object space to a Riemann surface, 
made of several copies, called sheets, of the Riemann sphere, cut and pasted. 
A point of the image space may then have several antecedents on the Riemann 
surface defining the object space. Example : two sheets for X >■ \J~X ^ Si countable 
infinity of sheets for x — > Log a;. 

As a consequence, to fill the reservoir of functions, one accepts all uniformiz- 
able applications, at the price of either restricting the object space by cuts, or 
defining a Riemann surface in the object space. 

Theorem. The general solution of a linear ODE is uniformizable. 

Proof. Let 

N A(k) 

E = Y,a k (x)^ r =0, a N (x) = l (3.1) 

fe=0 

be such an N th order ODE. Its general solution 

N 

u = y ^ Cjttj, Cj arbitrary constant, (3-2) 

has for only singularities those of the N independent particular solutions Uj, a 
subset of the singularities of the coefficients (|3(| §15.1). One knows where 
to make cuts or to paste the sheets of a Riemann surface in order to uniformize 
the general solution. QED. 

This has important consequences. Firstly, any linear ODE defines a func- 
tion (Airy, Bessel, Gauss, Legendre, Whittaker, ...). Secondly, in the needed 
reservoir of functions, one can put all the solutions of all the linear ODEs. 
Thirdly, a nonlinear ODE is considered as integrated if it is linearizable (of 
course via a finite linearizing expression). Fourthly, in order to extend the list 
of known functions by means of ODEs, it is necessary to consider nonlinear 
ODEs. 

Hence the problem stated by L. Fuchs and Poincare. 

Problem. Define new functions by means of ODEs, necessarily nonlinear. 

3.2 Second classification of singularities, differ- 
ent kinds of solutions 

There exist two features of nonlinear ODEs without counterpart in the linear 
case, they concern the location of singularities of solutions and the possible 
existence of solutions additional to the general solution. 

The singularities of the solutions of nonlinear ODEs may be located at a 
priori unknown locations, which depend on the constants of integration. 
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Definition (already given in the chapter Introduction). A singular point of a 
solution of an ODE is called movable (resp. fixed) if its location in the complex 
plane depends (resp. does not depend) on the integration constants. 

The point at oo is to be considered as fixed. A linear ODE has no mov- 
able singularities, the zeroes of its general solution depend on the integration 
constants and are sometimes called for this reason movable zeroes. 

Definition. The second classification of singularities is the distinction mov- 
able or fixed between singularities of solutions of ODEs. 

Among the four structures (critical or noncritical) and (fixed or movable) 
of singularities of solutions of ODEs, only one is an obstacle for this solution 
to be uniformizable and hence to define a function. This is the presence of 
singularities at the same time critical and movable. Indeed, in such a case, one 
knows neither where to make cuts nor where to paste the Riemann sheets, and 
uniformization is impossible. 

But let us come to the second distinction between linear and nonlinear 
ODEs : contrary to linear ODEs, nonlinear ODEs may have several kinds of 
solutions. 

Definition. The general solution (GS) of an ODE of order N is the set 
of all solutions mentioned in the existence theorem of Cauchy (section |5.3| ), 
i.e. determined by the initial value. It depends on N arbitrary independent 
constants. 

Definition. A particular solution is any solution obtained from the general 
solution by giving values to the arbitrary constants. A synonym in English is 
special solution. 

Definition. A singular solution is any solution which is not particular. Linear 
ODEs have no singular solution. 

Example. The Clairaut type equation 2it' 2 — xu' + u = has the general 
solution cx — 2c 2 , a particular solution x — 2, and the singular solution x 2 /8. 

A singular solution can only exist when the ODE 

E(u (N \u (N - 1 \... 1 u',u,x) = Q, (3.3) 

considered as an equation for the highest derivative u^ N ', possesses at least two 
determinations (branches), whose coincidence may define a singular solution. 
This is a generalization of the notion of envelope of a one-parameter family of 
curves. A practical criterium to detect the singular solutions will be given in 
section 5.1 . 

Painleve stated the following programme (|)0| p. 201, Oauvres vol. Ill p. 123; 
[|| p. 2, Oauvres vol. Ill p. 188) : 

"Determiner toutes les equations differentielles algebriques du premier ordre, 
puis du second ordre, puis du troisieme ordre, etc., dont Pintegrale generale est 
uniforme" . 

One notices that singular solutions are excluded from this statement. Indeed, 
they present no interest at all for the theory of integration, for, according to 
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the above theorem, they satisfy an ODE of a strictly lower order than the ODE 
under consideration and have therefore been encountered at a lower order in the 
systematic programme stated by Painleve. 

This problem (single valuedness of the general solution) splits into two suc- 
cessive problems whose methods of solution are completely different : absence 
of movable critical points, then absence of fixed critical points. Hence the final 
statement. 

Problem. Determine all the algebraic differential equations of first order, 
then second order, then third order, etc., whose general solution has no movable 
critical points. 

This class of equations is often denoted "with fixed critical points" . Let us 



prove that it coincides with the definition of the PP given section 1.4. Out of the 



four configurations of singularities (critical or noncritical) and (fixed or mov- 
able), only the configuration (critical and movable) prevents uniformizability : 
one does not know where to put the cut since the point is movable. 

We have now reached the usual definition, equivalent to the one of section 



1.4 



Definition. One calls Painleve property of an ODE the absence of movable 
critical singularities in its general solution. 



3.3 Groups of invariance of the PP 

In the fulfillment of the programme of Painleve, it is sufficient to take one 
representative equation by class of equivalence of the PP. There exist two re- 



lations of equivalence for the PP, defined in sections 3.3.1 and 3.3.2. Other 



relations of equivalence are defined in section 3.3.3, but they violate the PP. 



3.3.1 The homographic group 

Theorem. The only bijections (one to one mappings) of the Riemann sphere 
are the homographic transformations 

az + j3 

z — > -, aS — 7^ 0, (a, (3, 7, 5) arbitrary complex constants. (3.4) 

72: + S ' 

Proof. See any textbook. These transformations define a six-parameter 
group Tt called Mobius group, also denoted PSL(2,C). This group plays a fun- 
damental role in the present theory. Given two triplets of points, there exists a 
unique homographic transformation applying one triplet onto the other one. 

Theorem. The PP of an ODE E(u, x) — is invariant under an arbitrary 
homographic transformation of the dependent variable u and an arbitrary holo- 
morphic change of the independent variable x 

(u, x) ^{U,X):u= "fry + K X \ , X = t(x), aS-pj^O, (3.5) 
j(x)LJ + o[x) 
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where a, /3,7, S, £ denote arbitrary analytic (synonym : holomorphic) functions. 

Proof. Let xq be a regular point of (a, (3, 7, (5, £), and Xo = £(£0) its trans- 
form. In some neighborhood of xo, the transformation between u and U is close 
to a homographic transformation with constant coefficients and, according to 
the previous theorem, the first classification (critical, noncritical) is invariant : 
if Xq is critical (resp. noncritical) for u, then X is critical (resp. noncritical) for 
U, and vice versa. Since (a, (3, 7, S, £) do not depend on x$, the second classi- 
fication (fixed or movable) is also invariant. Thus the PP, which only depends 
on these two classifications, is invariant. QED. 

An element of this homographic group ( |3.5| ) will be denoted T(a, (3, 7, 5; £) 
or simply T(a,/3;£) in the case (7 = 0,<5 — 1). The representative equation is 
chosen so as to "simplify" some expression, e.g. a three-pole rational fraction 
the poles of which can be set at predefined locations like (00, 0, 1). 

Exercise 3.1 Choose a representative for the Riccati equation \l. \ ) in its class 
of equivalence under the homographic group. 

Solution. Two coefficients can be made numeric and the equation reduced to 

dU/dX + U 2 + S(X)/2 = 0, (3.6) 
under the linear transformation T(a!,/3;£) 

u = aU + (3, X = x, a = -l/a 2 , (3 = -(a' 2 + a ia2 ) / (2a 2 2 ) .D (3.7) 
This canonical form, in which S is called the Schwarzian, will be encountered 



again in section 6.3 



3.3.2 The birational group 

The PP is also invariant under a larger group (8^, [5^] , namely the group of 
birational transformations, in short the birational group, 

(u, x) (U,X) : u = r{x, U, dU/dX, d N ~ 1 U/dX N - 1 ) = 0, x = E(X), 

(3.8) 

(U, X) ^{u,x): U = R(X, u, du/dx, d^^^u/dx 1 ^- 1 ) = 0, X = £(x), 

(N order of the equation, r and R rational in U, u and their derivatives, analytic 
in x, X). 

For instance, given the ODE u" — 2u 3 = and the new dependent variable 
U = u' + u 2 , the algebraic elimination of (u 1 , u") among these two equations and 
the derivative of the second one yields the inverse transformation u = U' / {2U), 
which, once inserted in the direct transformation, yields the transformed equa- 
tion UU" - U' 2 /2 - 2U 3 = 0. 
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3.3.3 Groups of point transformations (Cartan equivalen- 
ce classes) 

The definition of to integrate as given page ^5] allows transformations outside 
the above two groups, which therefore may alter the PP. For instance, the 
unstable ODE 2uu' — 1 = is made stable by the change u 2 — ► U. One such 



group of point transformations, studied by Roger Liouville J82j, Tresse [109] and 
Cartan pi] , is defined as (it includes hodograph transformations) 

(u,x)+*(U,X) : u = f(X,U), x = g{X,U), U = F(x,u), X = G(x(&%) 

and the variables u and x are two equivalent geometrical coordinates. This geo- 
metric approach, in which provides a complementary insight to that of Painlcvc, 
which forbids to exchange the dependent and the independent variables, see sec- 
tion 



1.3 



The subgroup of fiber-preserving transformations 

{u,x)++(U,X) : u = f(X,U), x = g{X), U = F(x,u), X = G(#,10) 

whose equivalence classes are called Cartan equivalence classes, has been exten- 
sively studied by Kamran et ai, see e.g. ||| . 

3.4 The double interest of differential equations 

Let us return to the above problem. At each differential order of the pro- 
gramme, the results are twofold (this is the "double interest" of differential 
equations) : 

1. some new functions (defined from the general solution of a stable ODE 
which is not reducible to a lower order nor to a linear equation), 

2. an exhaustive list (i. e. a classification) of stable ODEs, which includes 
the ones defining new functions. 

Of course, each equation is characterized by one representative in its equiv- 
alence class. Thus, as seen in the introduction, the ODE for x(t) in the case 
b = 0, er = 1/3 of the Lorenz model is not distinct, under the homographic 
group, from the (P3) equation in the case a = j3 = 0, 7 = 5=1. 

For instance, the result for order one and degree one (the degree of an al- 
gebraic ODE is the polynomial degree in the highest derivative) is : no new 



function, one and only one stable equation which is the Riccati equation ([LI]) 

3.5 The question of irreducibility 
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The classical definition of irreducibility as given by the "groupe de ratio- 
nalite" of Jules Drach (Drach, in Cteuvres vol. Ill p. 14, @ p. 246, pjj ]) had 
some weaknesses pointed out by Roger (not Joseph) Liouville in a passionating 
discussion with Painleve in the Comptes Rendus (see Oeuvres vol. III). This is 
only after further mathematical developments, namely the differential Galois 
theory, that a precise definition of irreducibility could be given by Umemura 
| JllO| , sec Okamoto, this volume. This definition shares many features with 
the algorithm of Risch and Norman in computer algebra (which decides if the 
primitive of a class of expressions, e.g. rational fractions, is inside or outside the 
class) . 

3.6 The double method of Painleve 

To solve his problem, Painleve split it into two parts (this is the "double 
methode" |9§ p. 11) : 

• [1] (a local study) construction of necessary conditions for stability, 

• [2] (a global study) proof of their sufficiency \ either by expressing the gen- 
eral solution as a finite expression of a finite number of elementary func- 
tions (solutions of linear equations, . . . ), or by proving the irreducibility 
of the general solution and its freedom from movable critical points. 

The methods pertaining to each part are different. If some necessary condi- 
tion is violated in the first part, one stops and proceeds to the next equations. 
If one has exhausted the construction of necessary conditions, or if one believes 
so (indeed, this process, although probably finite, is sometimes not bounded), 
one turns to the explicit proof of sufficiency, i.e. practically one tries to integrate 
(no irreducible equation has been discovered since 1906). 

For a good presentation of ideas, see the book of E. Hille |]3| . 

Remark. The reason why movable essential singularities create difficulties 
lies in the inexistence of methods to express conditions that they be noncritical 
(Legons pp. 519 sq.). 

3.7 The physicist's point of view 

The physicist is not interested in establishing a classification nor in finding 
new functions. Usually, some differential system, whether ordinary or partial, is 
imposed by physics, and the problem is to "integrate" it in some loose sense. By 
the way, this loose objective is certainly the main responsible for the numerous, 
of course divergent, interpretations of "integrable" to be found in the physicists' 
world. 

The best applicability of the present theory arises when one knows nothing 
or very little on the possible analytical results : first integral, conservation laws 
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of PDEs, particular solutions 
section 



, Then, the first part of the double method of 



3.6 happens to be a precious integrability detector. We have already seen 



a rough version of it : the meromorphy test of section 2.3, the final version of 
which will be the Painleve test section 3.6. 

The loose objective of the physicist implies performing the test to its end, 
even if at some point it fails and should be stopped. One thus gathers a lot 
of information in the form of necessary conditions for a piece of local single 
valuedness to exist. This partial integrability detector can be called the "partial 
Painleve test" and will be examplified in section 3.7. 

Examining each condition separately, i. e. independently of the others, or 
simultaneously, one may then find pieces of global information like a first integral 
or a particular closed form solution. 
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Chapter 4 

The classical results 
(L. Fuchs, Poincare, 
Painleve) 



The problem of determining all stable equations has been completely studied 
for several classes of ODEs, while others are still unfinished. We review here 
the main results achieved to date. 



4.1 ODEs of order one 



The completely studied class is (L. Fuchs, Poincare, Painleve) 

E = P(u',u,x) = 0, P polynomial in (it', it), analytic in x. (4-1) 

When the degree is one, i.e. for the class vl — R(u, x) with R rational in 
u and an alyt ic in x, one finds one and only one stable equation, the Riccati 
equation (1.1). Since it is linearizable, this case defines no new function. 

When the degree is greater than one, one finds one and o nly o ne new func- 
tion, the elliptic function p of Weierstrass, defined from eq. ( 1.12 ). The stable 
equations are : all the ODEs whose general solution has an algebraic dependence 
on the arbitrary constant, plus five binomial equations with constant coefficients 
(i. e. u' n = P m (u), (m,n) 6 TV, P m polynomial of degree m) . Historically found 
by Briot and Bouquet, these binomial equations have the following solution (see 
e.g. || Table 1 p. 73) 



u' n = (u-o) n+1 (u-6) n_1 , n > 2, 



'-(X ~ Xq) 



(4.2) 
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u — a) 2 (u — b)(u — c), 



u — a)(u — b)(u — c)(u — d), 
(u — a)(u — b)(u — c)Y 



Acosh[B(x-x )] + C (4.3) 
— = Ap{x-x ,g 2 ,gz) + 5(4.4) 



1 



u — a 
1 



= Ap'(x-x ,0,g3) + B (4.5) 



1 



u- af(u-bf(u- c) 2 , 

u — c a — c 



Ap 2 (x-x ,g 2 ,0)(4.6) 



u-af{u-bf{u-c) A , 



u — a 



Ap 6 (x-x ,0,g 3 ) + B 1 (4.7) 



in which a,b,c,d are complex constants and (A, B,C, g 2 , g 3 ) algebraic expres- 
sions of (a, b, c, d). 
Remarks. 



1. Equation (4.3) is a degeneracy of (4.4). The T transformation u — » a^+u -1 
(<2i zero of P2n) generates eleven other equations with m < 2n, among 
them the Weierstrass equation (1.12). 



2. If the Weierstrass equation had not been known, it would have been dis- 
covered at this order one of the systematic process of Painleve. 



4.2 ODEs of order two, degree one 

The study of the class 

u" = R(u' ,u,x), R rational in u' , algebraic in u, analytic in x (4-8) 



was started by Painleve |9C|, |92j, 94 and finished by his student Gambier |56j . 

This class provides six new functions, the functions of Painleve, defined by 
the ODEs 



(PI) u" 
(P2) u" 

(P3) u" 
(P4) u" 
(P5) u" 
(P6) u" 



6u + x 

2u 3 + xu + a 

u a u' au 2 



•7« + -, 
u 



U X X 

U ' 2 3 ,3 , A m „,2 , o(J2 .\_. , P 



2u 2 u 



1 

2u 



1 



1 



u- 1 
1 



(u - I f 



1 



u U — 1 u — X 
u(u — l)(u — x) 



1 



x 2 (x — l) 2 



1 

x x — 1 
x- 1 



u p m(m + 1) 
7 - + 5^ -L, 

x u — 1 



1 



u — a; 
x(x- 1) 



(u — l) 2 (u — a;) 2 
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depending on respectively 0,1,2,2,3,4 complex parameters, since the nomogra- 
phic group allows to restrict to 7(7 — 1) = 0, 5(5 — 1) = for (P3), and to 
5(5 - 1) = for (P5). 



The stable equations (4.8) define 53 equivalence classes under the homo- 
graphic group, including of course the six above ones. They split into 50 with 
R rational in u and 3 with R algebraic in it, and their list can be found in : 
the original articles of Gambier @, ||, @, Gambier These 1910, Ince 1926 
[ |66| (caution : the numbers 5, 6, 48, 49, 50 of Gambier are changed to 6, 5, 49, 
50, 48 in Ince), Murphy 1960 fg, Davis 1961 @, Bureau M. I 1964, Cosgrove 
1993 [Q. In fact, the historical list of Gambier mixes two notions on purpose, 
namely the irreducibility and the homographic group, which makes this number 
53 rather arbitrary; for instance, the classes numbered 1 — 4, 7 — 9 by Gambier 
have been united by Gamier |38) into the single class 

u" = 5(2u 3 + xu) + 7(6it 2 + x) + f3u + a, (4.9) 

a stable equation admitting a second order Lax pair. 



The 50 stable equations (4.8) with R rational in u define 24 equivalence 



classes [p6| under the birational group and less than 24 Cartan equivalence 



classes |65j under the group of point transformations ( 3.10 ). 

This extremely important result (the discovery of six new functions, nowa- 
days frequently encountered in physics, and the exhaustive list of 50 + 3 equa- 
tions) deserves several comments depending on the field of interest. 

• (Practical usage). Given an algebraic second order, first degree ODE, 
either it is transformable by a T transformation ([T^) into one of the 
50 + 3 equations or not. If it is not, it has movable critical points. If 
it is, it is explicitly integrated and, by looking in the table of Gambier, 
its general solution is a known finite single valued expression made of the 
following functions : solutions of linear equations of order at most four, 
Weierstrass, (PI) to (P6). 

• (Movable singularities). The only movable (and of course noncritical) 
singularities of these ODEs are : poles for 50 + 2 of them, in addition a 
nonisolated essential singularity for 1 of them. This result of Painleve and 
Gambier (poles are the only movable singularities of stable second order, 
first degree equations rational in (it', it) and analytic in x) is often believed 
to be more general, leading to wrong definitions for the PP; it is no more 
true for third order, degree one, or even second order, degree higher than 
one. 

• (Fixed critical singularities). (PI), (P2), (P4) have none, (P3) and (P5) 
have two transcendental critical points (00, 0), both removable by the uni- 
formizing transformation x — > e x . (P6) has three transcendental critical 
points (00, 0, 1). 
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• (Dependence on the arbitrary constants). What characterizes the 6 Pain- 
leve equations among the 24 is the transcendental (i.e. not algebraic) de- 
pendence of their general solution on both constants of integration. The 
24 — 6 equations whose general solution does not involve (P1)-(P6) have 
either an algebraic dependence on both constants or a semi-transcendental 
dependence (algebraic for one, transcendental for the other one). 

• (Confluence). By a confluence process (Painleve fl94|] , Gambier These, see 
Mahoux, this volume), (P6) generates the five others and (PI) generates 
the Weierstrass equation, so up to now algebraic equations have only de- 
fined one master function. 

• (Monodromy). (P6) was found independently by R. Fuchs []53| and Schle- 
singer |l05|| in the twenty-first problem of Riemann. Given the second 
order linear ODE y" (t) 4- ai(t, x, u)y'(t) + a 2 (t, x, u)y — with four Fuch- 



sian singular points t = (oo,0, X,x) (see definitions section 5.2.1) and an 
apparently singular point u, the necessary and sufficient condition for the 
group of monodromy (see Mahoux, this volume, for definitions) to be in- 
dependent of x is that u(x) satisfies equation (P6). 



3 ODEs of higher order or degree 

Painleve's opinion was that no new function should be expected at third 
ler and that one should go to fourth order. In fact, despite huge efforts, no 
w function has yet been found. 

• ODEs of order two, degree higher then one. 

Only some subclasses have been studied, and their classification is nearly 
finished. See Chazy (These), Bureau (M. Ill), Cosgrove (1993 @|^,||). 

Those of degree two have the necessary form 
[u" + E u' 2 + E lU ' + E 2 } 2 

= F u /4 + iqu' 3 + F 2 u' 2 + F 3 u' + F 4 , (4.10) 

with (Ek^Fk) rational in u and analytic in x. Its binomial subset (E = 
Ei = E 2 = 0) is classified in Ref. @. 

The binomial subset (u") n = F(u',u,x) of equations of degree n > 3 is 
classified in Ref. |o). 

• ODEs of order three, degree one. 

The classification is nearly finished. See Gamier (These), Chazy (These), 
Bureau (M. II, Cosgrove 0. 
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• ODEs of order four, degree one. 

The classification is just started. See Chazy (These), Bureau (M. II). 

For an account of similar work on PDEs, see Ref . [B4| . 
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Chapter 5 

Construction of necessary 
conditions. The theory 



The reader only interested in using the Painleve test may skip this chapter, 
whose relevant parts will anyway be referred to in next chapters. By so doing, 
however, his/her confidence in the Painleve test will falter at the first encounter 
of one of the innumerable so-called exceptions, counterexamples, and so on, 
which are published every year. 

This chapter describes all the methods to build necessary conditions (NC) 
for the absence of movable critical points in the general solution. Most methods 
are analytic, and we unify their presentation by describing each of them as a 
perturbation in a small complex parameter e, to which can then be applied the 
theorem of perturbations of Poincare, itself a generalization of the existence the- 
orem of Cauchy. One of them is arithmetic and leads to diophantine conditions 
on the Fuchs indices of a linear differential equation. 

What we try to emphasize is the quite small amount of nonlinear features 
in these methods. Indeed, most of the information is obtained by well known 
theories concerning linear equations, whether differential or algebraic. 



The detection of singular solutions is first explained in section 5.1. The 



linear ODEs are then reviewed from the point of view of interest to us. Then 
we state the fundamental theorems at the origin of all methods. For comparison 
purposes, two equations are defined which will be later processed by all methods. 
Finally, we describe each method and apply it to the two examples. 

Unless otherwise stated, the class of DEs considered is made of DEs (2.30) 
polynomial in u and its derivatives, analytic in x, with (E, u) multidimensional. 



5.1 Removal of singular solutions 
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Since the PP excludes the consideration of singular solutions, one must dis- 
card them as early as possible. 

Let us give a practical criterium to detect singular solutions. 
Theorem. A necessary condition for a solution of an ODE to be singular is 
the existence of a common finite root to E = and its partial derivative 
with respect to u^ N K If E(u,x) = depends polynomially on the two highest 
derivatives after factorization of this polynomial existence condi- 

tion in w*^ -1 ) (called discriminant), it is necessary that the vanishing factor 
has an odd multiplicity. 

Proof. See e.g. Chazy (The se). The condition is not sufficient, and details 



and examples can be found in [104] chap. 10. 

Hence the method : compute the discriminant, factorize it, discard the even 
factors, test each odd factor to check if it defines a solution to the equation. 

Chazy (These p. 358) was the first to notice the absence of correlation be- 
tween the structure of singularities of the GS and of the SS. Here are such 
examples. 

Single valued GS, SS with a movable critical point (Chazy, These p. 360) 

(u'" - 2u'u") 2 + Au" 2 {u" - u 12 - 1) = 0, 
discriminant = — 16u" 2 (u" — v! — 1), 

c 2 - 4 

GS : u = e ClX+C2 /ci + x + c 3 , SS : u = C 2 - Logcos(x - CM5.1) 

4ci 

GS with a movable critical point, single valued SS (Valiron tome II §148) 

27uu /3 - I2xu' + 8u = 0, 

discriminant = -12 3 x 27 2 u 2 (27w 3 - 4a; 3 ), 

GS : u 3 = c(x - c) 2 , SS : u 3 = (4/27)x 3 . (5.2) 

Single valued GS and single valued SS (Painleve BSMF p. 239) 



ii 



" 2 + 4u' 3 + 2(xu' - u) = 0, 



discriminant = — 8(2u' 3 + xv! — u), 

GS : u = (l/2)i/ 2 - 2v 3 - xv, v" = 6v 2 + x, SS : u = Cx + 2C 3 .(5.3) 



5.2 Linear equations near a singularity 

Our only interest here is to decide about the local single valuedness near a 
singularity x = xo, put for convenience at the origin by a homographic trans- 
formation (x — > x — x or x — > 1/x according as xq is at a finite distance or 
not). 

These results are detailed in the course of Reignier, this volume. 



40 



Consider the most general linear system, put in a form solved for all first 
order derivatives (the canonical form of Cauchy) 









-( 


5 






x^- = AU, (5.4) 
dx 

with U a column vector of N components and A a square matri x r ational 
in x. This can be the representation of the general scalar ODE ( |3.l| ), with 
U k = x k u( k \k = 0,...,N- 1 and b k = x N ~ k a k , e. g. for TV = 3 



(5.5) 



Definition. The point x = is called Fuchsian iff all solutions of ( |5.4D have 
a polynomial growth near it. It is called nonFuchsian if at least one solution 
has a nonpolynomial growth. 

Example. For the ODE v! + ax n u = 0, n integer, a nonzero, the Fuchsian 
case is n > —1 and the nonFuchsian one is n < —2. The solution u(n) for 
n = —2,-1,0 is u(— 2) = e a l x ,u{— 1) = x~ a ,u(0) = e~ ax and its singularity at 
x = is respectively : an isolated noncritical essential point, a critical point or 
pole or zero (depending on a), a regular point. 

Remarks. 

1. This definition is the one of modern authors 0. It involves a property 
of the solutions, not of the coefficients of the equation. Fuchsian denotes 
at the same time a case with the solutions u — (x,x 2 ) (classically called 
regular point) and a case with u = (x , x ) (classically called singular 
regular point). The motivation for such a definition is the difficulty to 



recognize it on the matricial notation. While in the scalar case (3.1) the 
canonical form defined by setting a^-i = provides an easy criterium 
to decide about the nature of the singularity, in the matricial case the 
example 

'nix 0\ TT _ _(x»\ TT /(T 



a/*=(j^ oj< ui = v ;j' u2 = v ij ^ 

shows the difficulty to do so. 

2. We avoid the usual synonyms regular singularity and irregular singularity 
for their built-in conflict. 

Definition. Given a point xq, a fundamental set of solutions of a linear 
ODE of order iV is any set of N linearly independent solutions defined in a 
neighborhood of xq. 
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5.2.1 Linear equations near a Fuchsian singularity 

Definition. Given a Fuchsian point x — 0, the eigenvalues i of A(0) are 
called Fuchs indices. The indicial equation is the characteristic equation of the 
linear operator A(0) 

lim det (A(x) - i) = 0. (5.7) 
Near a Fuchsian point x = 0, there exist a fundamental set of solutions 



x x 



5^ tf (z)(LogaO J ', i = l,N (5.8) 

3=0 



in which the A,'s are complex numbers (the Fuchs indices), m, positive integers 
(their multiplicity), cpij converging Laurent series of x with finite principal parts. 
Series ( |5.8[ ) are the simplest examples of -0— series. 

The necessary and sufficient condition of local single valuedness of the general 
solution of the linear equation is : Aj all integer, no Log terms. 
In the scalar case ([O]), the indicial equation is 

= lim x N ~ i E(x, x l ) = lim det (A(x) - i) (5.9) 

x^O x^O 

= 6 (0) + h(Q)i + &2(0)*(* - 1) + . . . + b N {0)i{i - 1) . . . (* - N + 1). 



Theorem. Given a Fuchsian point of the scalar ODE fl3.l| ), necessary and 
sufficient conditions for the general solution to be locally single valued near it 
are 

o the N indices are distinct integers, 

o N(N — l)/2 conditions for the absence of logarithms are satisfied. 

Proof. See Hille book p3| or any other textbook. 

In the matricial case (5.4), these conditions are replaced by 

o the N indices are integers, 

o the multiplicity of each index i is equal to the dimension of the kernel of 
A(0)-i, 

o all conditions for the absence of logarithms are satisfied. 

The search for the no-log conditions can be achieved in one loop, by requiring 
the existence of a Laurent series extending from the lowest Fuchs index ii to 

+oo 

x ii J2 u j x3 ( 5 - 10 ) 

3=0 
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and containing N arbitrary independent coefficients; this is a finite process, 
which terminates when j + i\ reaches the highest Fuchs index. Consider for 
instance the third order ODE admitting the three solutions 

Mi = x~ 2 , u 2 = x~ 3 + ax~ 2 Logx, u 3 = x~ 4 + bx~ 3 Logx + ^-x~ 2 (Logx) 2 , 

namely 



u 


Ui 


U 2 


u 3 


u' 
















u" 


u'l 


< 


u'i 










v!" 


< 


< 


«3" 



= 0. 



Its three Fuchs indices —4, —3, —2 are simple, it is sufficient that j runs from 
to 2 with i\ = —4, the condition b = is found at j = 1 and the condition a = 
at j = 2. 

In some of the next sections, we will encounter inhomogeneous ODEs in 
which the rhs is itself a Laurent series with a finite principal part, so we will 
have in addition to express the single valuedness of a particular solution as well. 
This can be incorporated in the single loop described above, provided it starts 
from the smallest of the two values i\ and the singularity order of the particular 
solution, imposed by the rhs. 

5.2.2 Linear equations near a nonFuchsian singularity 

Near a nonFuchsian singular point x = 0, there exist N linearly independent 
solutions 

rrti 

eW'tx'^VijiziKLogxy, z i = x 1 '«, i = l,N (5.11) 
j=o 

in which the q^s are positive integers, Qi polynomials, Sj complex numbers 
called Thome indices, tpij formal Laurent series with a finite principal part. 
The question of local single valuedness of the general solution cannot be settled 
so easily, because formal series are generically divergent. 



5.3 The two fundamental theorems 

Theorem I (Cauchy, Picard). Consider an ODE of order TV, of degree one in 
the highest derivative, defined in the canonical form 

^ = Kb, u],ieC,ue C N . (5.12) 
ax 

Let (x , uo) be a point in C x C N and D be a domain containing (x , u ). If K 
is holomorphic in D, 
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• there exists a solution u satisfying the initial condition u(xo) = Uo, 

• it is unique, 

• it is holomorphic in a domain containing (xo,Uq). 

Proof. See any textbook. For delicate points on this classical theorem, see 
Legons p. 394. The contribution of Picard is to have moved the holomorphy 
property from the hypothesis to the conclusion. 

There exists an important complement to the theorem of Cauchy, due to 
Poincare : the Cauchy solution is also holomorphic in the Cauchy data. 

Remark. More practically, the canonical form can also be defined as 

A*.. 

= K[x,u,u',...,u (N - 1) }. (5.13) 



dx N 



The theorem says nothing whenever the holomorphy of K is violated, as in 
the following two cases. 

Case 1. du/dx = u/{u — 1), at uq = 1, a point of meromorphy for K. 
The only way to possibly remove this singularity without altering the structure 
of singularities is to perform a T transformation ( |3.5[ ). The homography T : 
l/{u—l) = V yields a new K, defined by dU/dx = —U 2 — U 3 , which is indeed 
holomorphic in C x C, now making the theorem applicable. In order to shorten 
the exposition, this step of an homographic transformation will be omitted in 
the whole chapters |E] and [| a nd o nly reminded for the synthesis of all methods 



into the Painleve test section 6.6 



Case 2. du/dx — \J 4(u — e\){u — e%j{u — e^), at uq — ej,j — 1, 2, 3, critical 
points for K. 

Example, du/dx + u 2 — 0, with the datum u = uq at x — xq. The Cauchy 



solution is represented by the (infinite) Taylor series u — uq ^2^^[—{x—xo)uoY, 
a geometric series whose sum depends on one, not two, arbitrary constants, the 
arbitrary location x\ = x^ — u^ 1 of the movable simple pole; it only exists locally, 
inside the disk of convergence centered at xq with radius | tio | 1 - This sum is 
also represented by the Laurent series [x — xi) _1 . One notices the enormous 
advantage of the Laurent series : it reduces to one term, and it has a much 
larger domain of definition (the whole complex plane but one point). 

Lemma (Poincare, Mecanique celeste p9| ) . Consider an ODE of order N, of 
degree one in the highest derivative, depending on a small complex parameter 
e, defined in the canonical form 

^ = Kb, u, el, xeC, ueC N , e e C. (5.14) 
dx 

Let (xq, u , 0) be a point in C x C N x C and D be a domain containing (xo, u , 0). 
If K is holomorphic in D, the Cauchy solution exists, is unique and holomorphic 
in a domain containing (xq, Uq, 0). 
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Proof. See any textbook. Note that K may be independent of e. 

Definition. Given x, the application u — > E(x, u) and some point u^°\ one 
calls differential of E at point u^ - 1 the linear application, denoted E'(x, u^), 
defined by 

Vv : E'( Xj u(°>)v = Ii» E(x,u(°> + Av) - E(x,u(°>) 

This notion is known under various names : Gateaux derivative, linearized 
application, tangent map, Jacobian matrix, and sometimes Frechet derivative. 
Definition. Given a DE E(x,u) = and a point uo, the linear DE 

E'(x,u'°»)v = (5.16) 

in the unknown v is called the linearized equation in the neighborhood of u^ -* 
associated to the equation E(x, u) = 0. 

This is precisely the equation auxiliaire ( 2.20| ) of Darboux. The auxiliary 
equation of a linear equation is the linear equation itself. 

Let us define the formal Taylor expansions 

+ 00 +OG 

u = e" u(n) , K = e " K(n) - ( 5 - 17 ) 

n=0 n=0 



The single equation (5.14) is equivalent to the infinite sequence 
dn(°) 

n = 0: — — = K(°) =K[a ) u(°' ) 0] (5.18) 
ax 

rln( n ) 

n>l: — — = KW =K'[a;,u( o ),0]u( n ) + R^fe u<°\ . . . , u ( "-%.19) 
ax 

At order zero, the equation is nonlinear. 

At order one, the equation, in the particular important case when K is 
independent of e, is the linearized equation (without rhs, RW = 0) canonically 
associated to the nonlinear equation. 

At higher orders, this is the same linearized equation with different rhs 
arising from the previously computed terms, and only a particular solution is 
needed to integrate. 

Theorem II (Poincare 1890, Painleve BSMF 1900 p. 208, Bureau 1939, M I). 
Take the assumptions of previous lemma. If the general solution of ( 5.14 ) is 
single valued in D except maybe at e = 0, then 

• e — is no exception, i.e. the general solution is also single valued there, 

• every xv- n > is single valued. 
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Proof. See BSMF p. 208. The main difficulty is to prove the convergence of 
the series. This theorem remains valid if one replaces "single valued" (Painleve 
version) by "periodic" (Poincare version) or "free from movable critical points" 
(Bureau version). 

Remarks. 



This feature (one nonlinear equation (5.18), one linear equation (5.19) 
with different rhs) is a direct consequence of perturbation theory, it is 
common to all methods aimed at building necessary stability conditions. 



The equations may be differential like ( 5.18 )-( 5T9| ), or simply algebraic. 



Moreover, all the methods which we are about to describe (except the one 
of Painleve) will reduce the differential problems to algebraic problems 
keeping the same feature, and the overall difficulty will be to solve one 
nonlinear algebraic equation, then one linear algebraic equation with a 
countable number (practically, a finite number) of rhs. 

The two theorems and the lemma express a local property, not a global 
one, therefore they cannot serve to prove integrability as defined page p5| . 
Conversely, they can be used to disprove the PP. In the same spirit, it is 
generally useless to try and sum the Taylor or Laurent series which will be 
defined. Indeed, these series only serve as generators of necessary stability 
conditions. The proof of sufficiency is achieved by completely different 
methods. 

The two theorems only apply to ODEs written in the canonical form of 
Cauchy. 



As a summary, the equations successively involved are 



the original DE E(x, u) = also called unperturbed DE because e will be 
introduced into the equation from the outside, 

the perturbed DE E(X, U,e), obtained from the preceding one by some 
transformation (x. u, E) — > (X, U, E, e) called perturbation, 

a canonical form (it is not unique) dXJ/dX = K(X, U, e) of the perturbed 
equation, also called abusively perturbed equation, 



the infinite sequence (5.18)-(5.19|) of equations independent of e. 



The methods described in next sections for establishing necessary stability 
conditions consist of building one or two perturbed equations from the original 
unperturbed equation, then of applying the theorem II at a point xq which is 
movable. This movable point can be either regular (method of Painleve) or sin- 
gular noncritical (all the others), which will require its previous transformation 
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to a regular point (by a transformation close to u — > u^ 1 ) for theorem II to 
apply. One is thus led to the equations (5.19), i.e. to one linear DE with a 
sequence of rhs. In order to avoid movable critical points in the original equa- 
tion, one requires single valuedness in a neighborhood of xq for : the general 
solution of the linear homogeneous equation, a particular solution of each of the 
successive linear inhomogeneous equations. 

One must therefore express that a very special class of linear inhomoge- 
neous DEs has a general solution single valued in a neighborhood of x . Their 
lhs (homogeneous part) is the linearized equation (equation auxiliaire) of a 
nonlinear equation which has already passed the requirement n = of the- 
orem II. The rhs (inhomogeneous part) of equation n depends rationally on 
{u( k \ k = 0, . . . ,n— 1} and their derivatives, all single valued near xo since the 
necessary conditions have been fulfilled until n — 1. 

For the coefficients of the homogeneous linear DE, the point xq will appear 
to be either a point of holomorphy (method of Painleve) or a singular noncrit- 
ical point (the other methods). In the latter case, both situations (Fuchsian, 
nonFuchsian, see sections 5.2.1 and |5.2.2| ) will occur. 



5.3.1 Two examples : complete (PI), Chazy's class III 

Example 1 ("complete (PI)") (BSMF p. 224, @ §14.311 p. 329, @ p. 267, 

+ + eu 2 + fu + g = 0, (5.20) 



E 



d 2 u 



dx 2 ' dx 

with (c, e, /, g) analytic in x, and e nonzero. This equation arises in the system 



atic study of class (4.8) and has led to the discovery of (PI) . 

Under a transformation T(a,f3;£) (|3.5| ), equation ( 5.20 ) is form-invariant 
(@ P- 267) 



d 2 U 
dX* 



£" a! 
c-^-2— 
£' a 



d£7 ^ 
5 dX + V 



f + c— + 2e/3 - — 
a a 



+ [g + ff3 + cp + e/3 2 - (3"] a^C 2 = 0, a? ± 0. 



U_ 

(5.21) 



This allows to assign simple predefined values to as many coefficients as gauges 
in T, i.e. three. For any value of (c,e,f,g) it is possible to choose for the 
coefficients of dU/dX, U 2 , U the values 0,6,0, and this requires solving two 
quadratures and one linear algebraic equation for (a, /3; £) 



(Log a)' = - 
5 



2e 



t i2 



6 







2e 



(5.22) 



[/ + c(Log a)' + (Log a)" + (Log a)' 2 ]. (5.23) 



Consequently, in what follows, one always assumes c = 0, e = 6, / = in (5.20). 
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Example 2 (Chazy complete equation of class III p2|). 

- u xxx + -{2uu xx - 3m 2 ,) + aiu xx + ciuu x + c u x + d 3 u 3 + d 2 u 2 + diu + d = 0, 

(5-24) 

where (a, a^, Ci, di) are analyti c in x , and a nonzero. This one has led Chazy to 
the discovery of his equation ( |l.l7|) . 

Under a transformation T(a, (3\ £), equation ( 5.24 ) is form- invariant. For any 
value of (a, di, Ci, di) it is possible to choose the values 2, 0, for the coefficients 
of UUxx ,Uxx ,U 2 , and this requires solving the coupled ODEs for (a,/?;£) 
[notation A = a' /a] 

2A' - A 2 + 2a- 3 [{a 2 Cl + 18ad 3 )A + a 2 d 2 + 9a' d 3 - 3a 2 ai d 3 ] = 0(5.25) 

= Y> P = °~ 2 ( 6aA + 3a ' - aa i). ( 5 - 26 ) 

i.e. one Riccati equation followed by two quadratur es. C onsequently, in what 
follows, one always assumes a = 2, ai = 0, d 2 — in ( 5.24 ). 
None of these two examples has singular solutions. 



Exercise 5.1 Show the impossibility to cancel d 3 in (5.24) by choosing a,p,^. 



5.4 The method of pole-like expansions 



This is a reliable version of the meromorphy test given in sections 2.3 and 
Consider a movable singular point xq of either the general solution or a 



particular solution. Since u(xq) is not finite, the theorems of section 5.3 cannot 



be applied. It is nevertheless immediate to check that the perturbation 
x = x +eX, u={EXy^(eX) n u {n \x), E = {eXf ^{eX) n E^ (as), 

n=0 n=0 

(5-27) 

in which the key point is the dependence of vS n > on x, not X, generates equations 



£j( n ) = which only differ from the algebraic equations Ej = defined by (2 
by the replacement of xq by x. The identification is even complete if \ is defined 
by Xx = 1 instead of x — x ~ x o , see Remark page ^6|. 

Fortunately, the method we are about to describe has been made by Bureau 



(1939) |13J an application of theorems I and II, as will be seen in section 5.6 . 
This method of pole-like expansions is the most widely used in Painleve analysis. 
Initiated by Paul Hoyer |64| and Sophie Kowalevski |f72| , f73[ , it has been for- 
malised by Gambier |36|, revived by Ablowitz et at M who applied it to wide 
classes of physical equations, extended to partial differential equations (PDEs) 



by Weiss et al. [ 1 1 1 1 (WTC), with technical simplifications by Kruskal 



and Conte |28L E9j. Painleve himself never used "le precede connu de Madame 
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Kowalevski . . . dont le caractere necessaire n'etait pas etabli" (Acta pp. 10, 83, 



Cteuvres III pp. 196, 269), see section 2.6 



We now rephrase the steps and generated conditions of sections 2.3 and 2.4 



so as to adapt them to the new objective : the PP. The expansion is denoted 



u = ]Tu^+p, u ^0, X ' = l. 

3=0 



(5.28) 



First step. Determine all possible families (p, uo). The necessary condition 
on (p,u ) is 

• CO. For each family not describing a singular solution, all components of 
p are integer. 

It there exists no truly singular family (at least one component of p negative), 
the method stops without concluding. 
Remarks. 

• Some components of Uo can be zero, or even some components of u. 

• To avoid missing some fa mily, one sho uld firstly put the ODE under a 
canonical form of Cauchy ( 5.12 ) or ( 5.13 ), so as to enumerate all the points 
u which make inapplicable the existence theorem of Cauchy, secondly for 
each such point build a transformed ODE under an homography making 
the point regular for the Cauchy theorem, thirdly determine families of 



the transformed ODE as in section 2.3 



• The derivative of order k of x p does not behave like \ p k if P is positive 
and p — k negative. 

Second step. For each family, compute the indicial polynomial detP. None 



of the conditions CI, C2, C3 of section 2.4 is required for the existence of 



the Laurent series since we also accept particular solutions and only exclude 
singular solutions. 

Third step. Unchanged as compared to section 2.4. The condition C4 is 
unchanged. 



The resulting expansion (5.28) thus contains as many arbitrary coefficients 



Uj as the sum of the multiplicities of the distinct positive indices, in addition to 
the arbitrary location of xq, associated to the index — 1. 

For indices which are not positive integers, the method says no thing , not 
even that they should be integers, and in such a case the expansion ( 5.28|) only 
represents part of the general solution, without indication about some possible 
multivaluedness in the missing part. 

Remarks. 
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The semi- infinite Laurent expansion (5.28) for u about the singular point 
xq is equivalent to an expansion for w -1 about a regular point, expansion 
however different from the Taylor one. This is used in Bureau's method 



section 5.6 



We prefer the terms "pole-like singularity" to "pole singularity" , for the 
actual singularity of the general solution may not be a pole, as shown 



by the example of Chazy's equation ( 1.17 ), for which it is a movable 
noncritical essential singularity. 

Index —1 also corresponds to an arbitrary coefficient but, since the general 
solution cannot depend on more than N such arbitrary coefficients, some 
renormalisation occurs. In the example du/dx + u 2 = already consid- 



ered in section 3.1, the Cauchy solution near the regular point xo can be 
reexpanded 

o 

u= (-uoVix-xoY- 1 (5.29) 

J=— oo 

[if the example were not so simple, this would be a doubly infinite Laurent 
series] so as to exhibit an arbitrary coefficient at index j = i = — 1, the 
only index of this too simple ODE. Note the "pole-like" singularity Xo, 
which is in fact an apparent, inessential singularity, in this case a regular 
point! 



5.4.1 The two examples 

Example 1. "Complete (PI)" eq. ( 5. 20] ). Already handled in section 2.3. 
Example 2. Chazy's equation (1.17). 

First step. The dominant terms are among —u"',2uu" — 3m' 2 , d^u 3 , hence 
two possible families 



0, V) = (-1, -4) u = -6 E = v!" - 2uu h 



3u r 



(p,q) = (-2, 



d?,u, 



3«0 



E 



-2m" + 3u' 2 + d 3 u 3 . 



The second family only exists if rfa = 0. See section |5.9| for a direct proof that 
d-i = is a necessary stability condition. 

Second step. The indicial polynomial of the first family is 

X - {l - i] [-dl + 2u oX - 1 d 2 x - Q{u oX - r ) x d x + 2(u QX - 1 ) xx }x t ~ 1 
= -(i - - 2){i - 3) - 12(i - T)(i - 2) - 36(i - 1) + 2(-l)(-2) 
= _(i + 3)(i + 2)(i + l), (5.30) 

and the indices are —3, —2, —1; the algorithm stops here, due to the absence of 
positive integer indices. 

Provided d 3 — 0, the second family has the indices —1,0. 
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Third step (only for the second family provided d% — 0). At the index 0, the 
condition C4 Qo — is satisfied and the algorithm stops. 

Exercise 5.2 Find the families and indices of the following equations. 

u" -2 = 0, u = (x - a) (x - b) (5.31) 



uu" — 2u a = 0, u = a(x — xo) 1 
u" + 3uu' + u 3 = 0, u ' ' 



x — a x 



(5.32) 
(5.33) 



5.4.2 Nongeneric essential-like expansions 



Just like (5.2g), the expansion 



12 u ^ 

-3=0 



j+P 



u ^0, 



(5.34) 



valid outside a disk centered at xq i.e. in a neighborhood of the point oo, is 
locally single valued. From this downward Laurent series, one could conceive a 
"method of essential-like expansions" quite similar to the method of pole-like 
expansions, in order to generate necessary stability conditions, this time from 
the negative integer indices only. 

However, for most equations, this method is not applicable. For instance, 
with the example — v!" + 2uu" — 3m' 2 + d^u 3 — (a subset of (1.17)), none of 
the two expansions ( 5.34 ) with p = —1 or p = —2 exists, unless ^3 = 0, which 
is not a reason to conclude that ^3 must vanish. 

It only applies to the very restricted class of equations with constant coef- 
ficients invariant under a scaling law (x, u) — > (kx, k p u), having at least one 
pole- like family with a negative integer index other than —1. Even then, its 
failure to detect the movable logarithm in num erous e q uation s which have one 
makes it of very little use. Such equations are ( 5.12S ), ( 5.144 ) or the equation 



u'" — 7uu" + llu' 2 = whose single family p = — 1, Uq = —2 has only negative 



integer indices (—6,-1,-1). 



5.5 The a— method of Painleve 



Consider an ordinary differential equation ([2. 30 ), a regular point xq (i.e. a 
point of holomorphy of the function K when (2.3C) is written in the canonical 
form (5.12)), define a small nonzero complex parameter (which Painleve denoted 
a) and the perturbation 



a^0 



x 



aX, u = a p a " u(n) : E = " q 12 ' 



'E 



(n) _ 



= 0, (5.35) 



71=0 



n=0 
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where p is a sequence of constant integers to be chosen optimally (see example 
below), q another sequence of constant integers determined by p, then apply 
theorem II to the equation for u(X,a). 
At perturbation order zero : 



all the explicit dependence of coefficients on X is removed, i.e. all coeffi- 
cients of the equation are constant, 

for a suitable choice of p, there only survive a few terms. 

the equation is invariant under the scaling transformation 
(X,u(°),E(°)) -► (fcX,fc p u(°),A; q E( )) (physicists call such an equation 
scaled, or weighted). 

Definition. The simplified equation (equation simplifies) associated to a given 
perturbation ( |5.35 ) is the equation of order zero E(°)(:ro,u( )) = in the un- 



known u<W(X). 

The simplified equation admits the one-parameter solutions Uq\x — Xq) p 
where Uq '' are constants. Its above properties usually make it easy to study. 

Definition. The complete equation (equation complete), as opposed to the 
simplified equation, is the equation it self (|2.30 ). 



The value a = is forbidden in (5.35), but theorem II takes care of that 



The constants p and q must be integers, chosen so as to satisfy the holomorphy 
assumption in the small complex parameter of the above lemma. Moreover, 
since a linear ODE has no movable singularity and since all successive equations 
E(") =0,n> 1, are linear, the only way to have movable singularities, in order 
to test their singlevaluedness, is to select simplified equations which are truly 
nonlinear. 

The successive steps of the a— method and the generated necessary condi- 
tions for stability are (BSMF p. 209 §7 and footnote 1) 

First step. Find all sequences p of integers satisfying the holomorphy as- 



sumptions of Theorem II for the perturbation (5.35). Retain only those defining 
a truly nonlinear simplified equation. For each sequence p perform the next 
steps. 



Remark. If the ODE (2.30) has degree one and order N, the holomorphy 
assumptions of Theorem II require that the highest derivative contributes to 
the simplified equation. 

Second step. Find the general solution u^ ^ of the simplified equation. 

• CO. Require u' ) to be free from movable critical points. 



The general solution v of the auxiliary equation (5.16) of the simplified 
equation is then (BSMF p. 209 footnote 1) 



dfe— - — , dk arbitrary constants, (5.36) 

c fc 
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and, since u^ -* has no movable critical points, v has no movable critical points 
cither (theorem Legons p. 445). 

Third step. For each n > 1, define u'™-* as a particular solution of equation 
E(«) = o (linear with a rhs), by the classical method of the variation of constants. 

• CI. Require each u*-™- 1 to be free from movable critical points. 

These steps amount to require stability for all the sequence of perturbed 
equations, exactly as formulated in Theorem II. 
Remarks. 

• In the second step one can take for u^ ) either the general solution or a 
particular one, but not a singular solution. The drawback of a particular 
solution will be a lesser number of generated necessary stability conditions. 
This may be useful when the quadratures of third step are difficult with 
the general solution and easy with a particular solution. 

• At order n — 1 equation E^ 1 ' = may contain a rhs, making it different 
from the auxiliary equation. 

• If in the second step only a particular solution has been chosen, it is better 
that u^ 1 ) be taken as the general solution of equation E' 1 ' = 0. 

Many people have intuitively used the a— method, let us give two recent 
examples. 

Example 1. In the Lorenz model ( |l.3| ), the simultaneous change of variables 
(x, y, z) — > (£, 7], C) and parameters (6, a, r) — > (6, a, e) defined by 

S, — ex, rj — e 2 ay, ( = e 2 az, e 2 ar = 1 (5.37) 



led Robbins 102 to believe to have found a new integrable case, defined by 

cr^0,£ = 0,r = oo: first integrals £ 2 - 2C, + rj 2 + C, (5.38) 

while in fact the new dynamical system is just the simplified of the original one, 
integrable by elliptic functions. 

Example 2. The transformation t — * t 2 Log t with "t — ■> 0" from the Lorenz 
model to the system (24abc) of Ref. is in fact the a— transformation 
(x, y, z, t) = (e~ 1 X, s~ 2 Y, e~ 2 Z, t + eT) resulting in the system 

dX , , AY dZ , 

— =a(Y- X), — = -Y-XZ,—= XY, (5.39) 

whose general solution is elliptic. 
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5.5.1 The two examples 

Example 1. "Complete (PI)" eq. (|5~20l ) (BSMF p. 224 §15). The Cauchy 
form of the perturbed equation is 

aP -2 E _ - d ^ +--- j +6^+ 2 [^°) + . . .} 2 +a^[g(x Q )+aXg'(x ) + . . .] = 0. 

(5.40) 

First step. The holomorphy requirement — 2 < p < 2,p £ Z, selects five 
values for p, and the requirement of a nonlinear simplified equation only retains 
p = —2, i.e. q = p — 2 = 2p = —4. The g(xo) term, which could vanish, does 
not contribute to the simplified equation 

E (o) ^_^!i + 6u (o) 2 = . (5 . 41) 
dX z 

Second step. The general solution of this particular Weierstrass equation is 

u (0) = p(X - c , 0, g 3 ), (co, 53) arbitrary. (5.42) 

The auxiliary equation of the simplified equation is a Lame equation 

,2 

E'(X,u^)v = + 12p(X-co,0,g 3 )v = (5.43) 

whose general solution is a linear combination of dp/d g3 and dp/d co 

v = Cl (Xp' + 2p) + c 2 p', (5.44) 

without any movable critical singularity. 

Third step. The successive linear equations with their rhs are 

E'(X,u i0) )uW = (5.45) 
E'(X,u^)uW = -6u« 2 (5.46) 
E'(X,u^)u^ = -l^V (5.47) 

i?'(X,J )y 4 ) = -m^^^-G^ 2 -.go (5.48) 
(°)) U < 5 ) = -12( U « U ( 4 ) W 2 V 3 ))- 3o X (5.49) 

£'(xy°V 6) = -12(u( 1 ^ (5 )+ M ( 2 )u( 4 ))-6u( 3 ) 2 -i 5o 'X 2 (5.50) 
with the particular solutions 

u {n) = 0, n = 1,2,3 (5.51) 
«W = |[2Xpp' + 2p 2 -Cp'] (5.52) 

u (6) = || [2X 2 pp' + 2Xp 2 + (Xp' + 2p)C] (5.53) 
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,(6) 
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[(Xp' + 2p)(X 2 + 2X( - 2 Log a) 



+ {X 3 p + X 2 (-2XLoga + 2 J LogadX)p'], (5.54) 

where the functions £ and a obey £' = — p,<r' = C°\ To prevent movable 
logarithms at n = 6 it is necessary that g"(xo) = 0. Since xq is arbitrary, this 
condition is Vx : g"(x) = 0, and Painleve proved it to be sufficient, thus defining 
the (new in the sense of Section 3.2) function (PI) with the choice g = x. 
Remarks. 

• The reason why ttW , , can be chosen zero is given on page The 
reason given in ref. J75| p. 120 is not correct : even if the general solutions 
uW, u( 2 ), it( 3 ) are meromorphic, they can in principle (this does not occur 
for the ODE under study) generate some multivaluedness further up in the 
computation. The theorem proven in Lecons p. 445 is quite profound : if 
a [second order in Legons] ODE is stable, its general solution has a single 
valued dependence on the integration constants. 



• Taking the particular solution u^ ' 
one p (see first remark section ^7. 



l/(X — cq) 2 instead of the general 
makes all computations immediate 
(C = l/(X — Co), a = X — Co). For this particular equation, one would not 
miss the generation of the only necessary stability condition. 



Exam-pie 2 (Chazy complete equation of class III ( 5.24 )). 

First step. For the canonical form of Cauchy of the perturbed equation 



d 2 (u 



(o) 



AX 2 

-a 3 - p d (x ) + 



+ 2aP+ V°> 



= 0, 



dV°> 



dx 2 



- 3a p 



dV°) 



dX 2 



the holomorphy condition is — 1 < p < 3,p G Z, which the condition for a truly 
nonlinear simplified equation restricts to p = —l,q = —4. The value p = —2 
[B2l pi) of the method of pole-like expansions is therefore forbidden. 



Second step. The simplified equation is that of Chazy (1.17), whose general 
solution u' ' is ]l7|, [l8| an algebraic transform (finite single valued expression) 
of the Hermitc modular function y(X) 



u 



(0) _ 



[Log^zr 2 ^-!)- 2 )] 



(5.55) 



evaluated at the point (ciX + C2)/(c3X + C4),ciC4 — C2C3 = 1 and thus obviously 
depending on three arbitrary constants. 

The auxiliary equation of Chazy's simplified equation 



E^'v = [-d% + 2u^d 2 x - &ufd x + 2i 



(0) 1 

xx\ 







(5.56) 
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has the three independent solutions du^ /d Ci , z = 1,2,3, all single valued. 

Instead of the general solution u^°\ which would make the computations 
rather involved, let us restrict to the two-parameter particular solution 

= — Gx 1 + c X 2 i X = X — xq, (xo, c) arbitrary constants, (5.57) 

for which the auxiliary equation admits the general solution 

v = k 2X - 2 + k 3X ~ 3 + M4, Vi = c~ 2 X - 2 (e~ 2c ^ - 1 - 2c X 1 ), (5.58) 

with (&2, k 3 , hi) arbitrary constants. 

Third step. The successive linear equations with their rhs are 

-#'uW = ciW 0) w (0) ' + rf 3 W° )3 (5.59) 
-E^'u^ = 2u^u^" - 3uW 2 + c^Xu^uW + c lfi u^' + u^) 
+d 3A Xu (0)3 + 3d 3 , u (0)2 u (1) + c , u (0) '- (5.60) 

A particular solution of the first one is provided by the method of variation 
of the constants = K 2 {X) X ~ 2 + K 3 {X) X - 3 + K^X)^ 4 

K' 2X - 2 + K' 3X - 3 + K' iVi = (5.61) 
-2K' 2X - 3 - 3K' 3X - A + KM = (5.62) 
6K 2X - 4 + 12K' 3X - 5 + K' iV 'l = c 1 ,o(-6 X - 1 + cx- 2 )(6x- 2 -2cx- 3 ) 

+rf 3 ,o(-6x~ 1 + c X - 2 ) 3 . (5.63) 

To prevent a movable logarithm in K 2 (resp. K 3 ), it is necessary that, in the 
rhs of last equation, the coefficients of X ~ 5 and X ~ 6 vanish : 

V(z , c) : -2c 2 ci(x ) - 18c 3 rf 3 (x ) = 0, c 3 d 3 (x ) = 0, (5.64) 

hence the two necessary stability conditions Wx d 3 (x) = c\(x) = 0, obtained at 
the perturbation order n = 1. We leave it as an exercise to check that, after 
completion of n — 4, one has obtained all the conditions (c\ = cq = d 3 = d\ = 
do = 0) which Chazy proved to be necessary and sufficient. 

Theorem. For any family of the method of pole-like expansions, the value 
i = — 1 is a Fuchs index. 

Proof. Let u ~ u x p be such a family and E(x, u) be the dominant terms. 
The equation E(xq, u) = admits as a particular solution the monomial X — ► 
u = Uo(xq)(X — Xq) p , therefore the linearized equation at the leading term 
(2.2C) admits as a particular solution its derivative with respect to Xo : X — > 
const x dx (X — Xo) p . Since at least one component of p is negative, the 
associated component of dx (X~Xq) p is proportional to (X — Xo) p_1 , therefore 
i = — 1 is a root of the indicial equation ( 2.36| ). □ 
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5.5.2 General stability conditions (ODE of order m and 
degree 1) 

Using his method, Painleve could obtain quite general necessary stability 
conditions for algebraic ODEs of arbitrary order and degree, cf. BSMF p. 258, 
Acta p. 74, Chazy (These). Consider the class, defined in the canonical form of 
Cauchy 

U M =R{u { - m - 1 \u { - m - 2 \...,u',u,x), (5.65) 

with R rational in u and its derivatives, analytic in x [for R algebraic, and for 
arbitrary order and degree, cf. Acta pp. 73, 77]. Necessary stability conditions 
are : 

CI. As a rational fraction of u^ m ^ x \ R is a polynomial of degree at most two 

U M = Au {m ~ 1)2 + Bu {m - 1] + C. (5.66) 

C2. As a rational fraction of u^ m ~ 2 \ A has only simple poles a,i with residues 
Ti equal to 1 — 1/m, rn nonzero integers possibly infinite 



i - iM 

,(m-2) _ 



(5.67) 



The above sum is finite. 

For second order m — 2 the fraction A has at most four simple poles, and 



the set of their residues can only take the five values of Table 

4 4 



.4 



^ u — a; ' ^ ' 



1 

2, Ti = 1 , ni g Z or n% 

m 



oo. 



(5.68) 



Table 5.1: Order two, degree one. Number of poles (nonzero ri), list of their 
residues. The poles may be located at oo and may not be distinct. The type 
numbering convention is that of (H| Table I p. 169). The least common mul- 
tiplier (1cm) is shown for convenience. 



Type 


lcm(ri) 




ri 




7"4 


I 


n > 1 


l + l/n 


1 - 1/n 








III 


2 


1/2 


1/2 


1/2 


1/2 


IV 


3 


2/3 


2/3 


2/3 





V 


4 


3/4 


3/4 


1/2 





VI 


6 


5/6 


2/3 


1/2 






Note the one-to-one correspondence between Table 
powers of the five Briot-Bouquet equations page M. 



5.1 



and the list of 
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Exercise 5.3 For the six equations (Pn), determine the set (oj, r<) of 
simple poles with their residue. 



Solution. 



(P6) (oo,l/2), (0,1/2), (1,1/2), (a:, 1/2), (5.69) 

(P5) (oo,l/2), (0,1/2), (1,1), (5.70) 

{PA) (oo,3/2), (0,1/2), (5.71) 

(P3) (oo,l), (0,1), (5.72) 

(P2) (oo,2), (5.73) 

(PI) (oo,2). (5.74) 



For instance, (P4) belongs to type I of Table 5.1, and it is also a confluent 
case of types III, V, VI. □ 

The similar finite lists of admissible values of A for any order m can be 
found in Painlcvc (m = 3 Acta p. 68, Oeuvres vol. Ill p. 254; m > 4 Acta 
p. 75, Oeuvres vol. Ill p. 261). 

C3. As rational fractions of u^ m ~ 2 \ B and C have no other poles than those 
of A, and these poles are all simple. Writing B, C as rational fractions 
of it( m ~ 2 ) whose denominators are that of A, this implies the degrees 
limitations 



(order 2, degree 1) : deg num B < 1, deg num C < 3. (5.75) 

C4. (Chazy, These). Every ODE u (™- 2 ) - ai = (a denominator of A) is 
stable. 

C5. ([^i)). All polynomial degrees in u^ k \ k = 0, . . . , m — 2, (of the numerator 
and denominator of A, B, C written as irreducible fractions of u and its 
derivatives) are limited, except in the "Fuchsian" case n.i — —2,ri = 3/2 
(see Ref. §§ for details). 

For additional conditions, see Ref. [pij . 



5.6 The method of Bureau 

Firstly, this method exhibits a linear differential equation with a Fuchsian 
singularity which allows to interpret the indices i in the recursion relation of 
Kowalevski as Fuchs indices. Secondly, it brings rigor to the heuristic method of 
Kowalevski and Gambier. However, the generated no-log conditions are identical 
to those of the method of pole-like expansions. 
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Consider an N th order ODE E(x,u) = (for simplicity, one assumes u and 
E unidimensional; the multidimensional case is handled in [|l9|) and a movable 
noncritical singular point xq where the general solution behaves like u ~ uq(x — 
xq) p , with p a negative integer to be determined. 

The integer p is computed by the method of pole-like expansions and the 
highest derivative is required to contribute (M. II p. 9) in order to be sure that 
one deals with the general, not a singular, solution. 

One wants to apply the two fundamental theorems. Since the singularity 
xo violates the holomorphy assumption of theorem I, one defines an equivalent 
differential system (in fact two systems) for which xq is a point of holomorphy. 
These systems will depend on a perturbation parameter e. 

One first defines two new dependent variables (z, U) by the relations (Gam- 
bier These p. 50, Bureau 1939) 

dz 

u = sz p , — = 1 + Uz, s ^ 0. (5.76) 
da; 

Elimination of u and the derivatives of z (M. II pp. 13, 77) 

z~ p u = s (5.77) 

z-p+^P = pa+ (?± + p a u)z (5.78) 



dx \dx 

,2 



P+2 d^ = p(p-i)«+(2p^+P(2p-i)«tr) 



^l +2 ppU+p 2 sU 2 +ps^-)z 2 (5.79) 
dx^ dx da; / 

etc, transforms E into 

E^E( X ,U,-,...,^ r ,s,-,...,^ r , Z )=0, (5.80) 

an equation for U of order N — 1 polynomial in z. For the equivalent system 
(|7r|), ( ^8C| ) made of two ODEs of orders one and N — 1 in the unknowns 
(z, U), the point z = is still a point of meromorphy, see examples below. 

To remove it, one introduces a dependence in a small nonzero parameter e 
to obtain a perturbed system to which Theorem II can be applied. Two such 
perturbations have been defined (Bureau 1939). 

First perturbation of Bureau 

+ OG 

x = x +eX, z = eZ, U unchanged : E = (sZ) qS ^{eZ) n E^ = 0, (5.81) 

?i=0 

where the positive integer —q is the singularity order of E. The coefficients must 
be expanded as Taylor series like in the a— method 



s(x) = s + (eX)s' + ..., s [ k) = ^Oo), a(x) = a + {eX)a' + .... (5.82) 



d (k) 



dx 
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Expansions up to order one in e for the above derivatives ( |5.77 )-( 5.79| ) are 

z -p u = s + {s' X)e + 0(e 2 ) (5.83) 
z~ P+1 ^ = ps + ((ps' Q )X + (s' +p S0 U)Z)e + O(e 2 ) (5.84) 



+P Up - 1)4* + (H + (2p - l)s Q U + *oZ^)Z 



0{e 



2 > 



etc, together with dZ/dX = 1 + eZC/ = 1 + 0(e). 

Order zero is an algebraic equation E(°'(xq, so) — for the nonzero coeffi- 
cient So- 

Order one is subtle : it filters out all terms nonlinear in U and its derivatives 
d (k) U/dX k , and it extracts the contribution of d (k) U/dX k from the term z + 
in the expansions ( 5.77 )-( 5.79| ). This results in 

X ^ d (k) U 

E {1) = A— + B+Y,c k Z k = 0, (A, B,c k ) constant. (5.86) 

fe=0 

Since dZ/dX is unity at this order, this is a linear inhomogeneous ODE of order 
at most N — 1 for U with constant coefficients, whose homogeneous part is by 
construction of Fuchsian type (exactly one singular point Z = 0, of the singular 
regular type) and even Eulerian type. 

In order to be sure of dealing with the general solution of the original nonlin- 
ear ODE, the linear ODE ( 5.86 ) must have exactly the order N — 1; a necessary 
stability condition for the nonlinear ODE is the single valuedness of the general 
solution of the linear ODE (5.86). Hence the necessary conditions, for each 
value of (p, so) 

- the order of the linear ODE at perturbation order one is exactly N — 1 ; 

- its TV — 1 Fuchs indices are distinct integers; 

- if is an index, the rhs vanishes (A = B = condition for the particular 

solution to contain no logarithm). 

Since Q5.86 ) is Eulerian, these conditions are sufficient for the general solution 
of the linear ODE (5.86) to be single valued, but only necessary for the stability 
of the nonlinear ODE. 

Higher perturbation orders yield no information. The reasoning is then that 
any condition thus found at x = xq, such as sq = 1, is valid at any x since xq is 
arbitrary. 
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Second perturbation of Bureau 

+00 +OC 

x unchanged, z = eZ, U= ^{eZ^U^ : E = {eZ) q ^{eZ) n E^ n \ 

n—1 n— 

(5-87) 

Expansions for the above derivatives ( 5.77 )-( 5.79 ) are 



z r u 
dx 

z 2-pd_U 

da; 2 



ps + [psU {1) + s') eZ + psU (2) {eZf + P sU (3) (eZ) 3 + 0(e 4 ) 

p{p - l)s + p ((2p - l)sU {1) + 2s') eZ 

+ (p 2 s(2UW + [/(i) 2 ) + 2 ps'U^ + s" + P8^£p) {eZf 



+0(e 3 ) 



etc, together with edZ/dx = 1 + eZU^ + {eZfU^ + 0(e 3 ). 

Equation (x, s) = 0, s 7^ 0, is the same algebraic equation as above for 
the unknown s(x), not s(xq). Each perturbation order n > 1 dehnes a linear 
algebraic equation 

Vn > 1 : P{n)U {n) + Q n (x, C/ (1) , . . . , U {n - 1} ) = 0, (5.88) 



where P(n) is the indicial polynomial of Fuchsian equation (5.86), and Q n de- 
pends on the previously computed coefficients. Necessary stability conditions 
Qi = arise at every value of i which is also one of the N — 1 Fuchs indices. 
These conditions ar e iden tical to those of the method of pole-like expansions, 
as proven in section 



5.6.1 



The successive steps and generated necessary conditions of the method of 
Bureau are 

Step a. Determine all possible p like in the method of pole-like expansions 
(details M. I p. 256, M. II p. 9). For all p satisfying (CO, CI), perform 
step b. 

CO. All p are integers. 

CI. The linear ODE (order one of first perturbation) has exactly order N—1 
[this holomorphy condition excludes for instance p = —2 in Chazy]. This 
implies the necessity for the highest derivation order to contribute to the 
dominant part during the computation of p. 

Step b. Solve the algebraic equation for so at order zero of first perturbation. 
For all nonzero sq perform steps c and d. 
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Step c. Solve the linear inhomogeneous Euler equation for U(Z) at order one 
of first perturbation. 

C2. Its N — 1 Fuchs indices are distinct integers. 

C3. If is an index, the inhomogeneous part vanishes. 



Step d. Solve the linear algebraic equation ( 5.88| ) (order n of second perturba- 



tion) from n = 1 to the highest Fuchs index. 
C4. Whenever the order n in step d is a Fuchs index i, Qi is zero. 

As compared with the a— method, these stability conditions are directly 
taken at x, not at xq- However, the method provides no conditions from the 
negative integer indices. 

5.6.1 Bureau expansion vs. pole-like expansion 

Let us first prove the existence of a one-to-one correspondence between the 
coefficients U^> of Bureau (second perturbation) and those Uj of the method 
of pole-like expansions. The relations defining Bureau coefficients are 

u = sz p (5.89) 

^ = l + U^z + U^z 2 + 0(z 3 ) (5.90) 
dx 

and those defining the pole-like expansion are 

u = X P {u + u 1 x + U2X 2 +0( X 3 )) (5.91) 

^ = 1. (5.92) 
dx 

The property \x = 1 of x nrs t ensures s = uq [taking \ — x ~ xq would just 
create useless complications]. The elimination of u between ( 5. 89] ) and ( 5.91 ) 
yields 

x(i + ^x + ^x 2 + o(x 3 " 

\ U Uq 

- x(l + ^x+ 2pU2+ 9 { *- pH X< + 0( x <)). { r,<m 

\ pu 2p z uf ) 
Let us invert this Taylor series z of x m t° a Taylor series x °f z 

x = z(i-^z+ - 2pU2 ^\ +p)u * zi + 0(z*)). ,r,c,j, 



pu 2p 2 UQ 



02 



One finally substitutes this x an d dz/dx, both Taylor series in z, into eq. (5.92) 
to obtain 



dx 



= ( l - + 3 ~ 2PM2 + (3 + P)ul z* + 0(z 3 ) 



xfl + f/Wz + plV + Ofz 3 ))-— ( ^) z 2 + 0(z 3 ).(5.95) 

pdx \uqJ 

The identification of the lhs and rhs as series in z provides the correspondence 
between the two sets of coefficients 

s = uq (5.96) 

= *± (5.97) 
pu 



jjmq \pu Q J 2p z UQ pdx \u 



(5.98) 



or 



"n 



Ul 



u 2 



s 
2 

O 



p- 



ip + l 



(5.99) 
(5.100) 

(5.101) 



24 6" dx 

This bijection between the coefficients induces a bijection between the equa- 
tions #(«) = of the expansion of Bureau (second perturbation) and the equa- 
tions Ej = of the method of pole-like expansions, hence a bijection between 
the no-log conditions. 

This proves the equivalence between the method of pole-like expansions and 
the second perturbation of Bureau. As to the first perturbation of Bureau, it 
brings quite important information not obtaibable by the method of pole-like 
expansions. 



5.6.2 The two examples 



Example 1. "Complete (PI)" eq. ( gjCj ). See Bureau M. I eq. (17.3), (21.2). 
The unperturbed equivalent meromorphic system ( 5.76| )-( |5~80" ) in (z, U) is, 
in Cauchy form 

dz 

1 + Uz (5.102) 



da; 

^ = 6«V-(p-l) 
da; 



z- 2 - 



2- + {2p-l)U 
s 



s'pU 



ps 



ps 
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Step a. The only value is p = —2, integer. The original ODE then reads, by 
increasing powers of z 



E 



6s(l 
-2s 



2s 



2 ds 
s dec 
1 d 2 s 



dU 2ds u .i. ' l 2 
dx s da; 2s dx 2 



z- 2 - g(x) = 0(5.103) 



Step b. Equation E^°> = 6so(l — so) = has for only nonzero solution so = 1. 
Step c. At order one 

f£ _ + K, - 2 4 + ,4§ _ 0. g - 1 + <**). (5.104) 

2So QA So So Z OA 

The only Fuchs index is i — 5. There is no condition on the rhs. 
Step d. The computation presents no difficulty. 



1, [/(l) = [/W = [/O) = 0, 



I 7-/(5) 

4' 



<?'(*) 



Qe 



(5.105) 



Remark. On the Cauchy form ( ^>. 102 ) with p = — 2, s = 1, one sees easily 
how perturbations I and II remove the meromorphy. 

Example 2 (Chazy complete equation of class III ( |5.24| )). 

Step a. The two solutions are p = — l,p = —2. For p — —2 the computation 
of the linear equation (5.86) yields a zero coefficient for d 2 U/dX 2 , thus 
violating condition CI. 

For p = — 1 the original ODE then reads, by increasing powers of z 



E 



s(s + 6)z~ 



„ rT 6 ds „ ds , , 

6Z7 — + 2— - Cl s + d 3 s 2 

s dx dx 



„/ ^dfj 9.ds .,_ 

2 a + 2)— + ((2 - - — - c lS )U - c 
dx s dx 



^ ds 3 ds ^ 2 ^ ,^ _|_ 3^ d 2 s 
dx s dx s dx 2 



+s 



d 2 C/ 



dx 2 



,dsdC7 ,3ds 
dx dx 



m.. <„ •. .. rr Co ds 

3 ^ + ^" c ^ + dl + 7d^ 



1 d 3 s 3 ds . 



. .. , dU , 
-— - -~rU 2 - 3U— + C/ 3 

s dx J s dx dx 



+ d . (5.106) 
Step b. Equation E^ = Sq(sq + 6) = has for only nonzero solution sq = —6. 
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Step c. At order one 



^d 2 U _d?7 

s a 3i0 



Z 2 -r^ - 2(s + 2)Z— + 1217 - s ci,o - sgd 3 

So OA ^ dA 



+ (2 - -)s + (2 + -)s ^ = 0, ^ = 1 + 0(^.107) 
so so Z dA 

The Fuchs indices are i — —4, —3, there is no condition on the rhs, and 
the algorithm stops here, due to the absence of positive integer indices. 

5.7 The Fuchsian perturbative method 

It allows to extract the information contained in the negative indices |fj2f , 
thus building infinitely many necessary conditions for the absence of movable 
critical singularities of the logarithmic type |jl| . 

The perturbation which describes it is close to the identity 

+ oo +oo 

x unchanged, u = ^ e n u (n) : E = ^ e"E (n) = 0, (5.108) 

n=0 n=0 

where, like for the a— method, the small parameter e is not in the original 
equation. 



Then, the single equation (2.30) is equivalent to the infinite sequence 



n = 0:E (0) ee E(x,u (o) ) = (5.109) 
Vn>l:E< n > = E'(x,u^)u^+R( I1 )(x,u^,...,u( n - 1 )) = 0(5.110) 

with R/ 1 ) identically zero. From Theorem II, necessary stability conditions are 
- the general solution xv-°> of (|5.109|) has no movable critical points, 



- the general solution of ( 5.110] ) has no movable critical points, 



- for every n > 2 there exists a particular solution of ( [5. 110 ) without movable 
critical points. 

Order zero is just the complete equation for the unknown u^ ) , so, to get 
s ome i nformation, one must apply Theorem II for a perturbation different from 
( [5.108 ). Since Bureau has proven that the method of pole- like expansions, with 



more rigorous assumptions, can be casted into an application of the two basic 
theorems, one uses it at order zero, only to obtain the leading term u' ' ~ Uq x p 
of all the families of movable singularities. 

First step. Determine all possible families (p,^ "*) 

u<°> ~ u V, E<°) ~ E< V, u[°» + (5.111) 
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which do not describe a singular solution, by solving the algebraic equation 

E< 0) = limy-^ukV) = 0. (5.112) 

• CO. All components of p are integer. 

If there exists no family which is truly singular (at least one component of 
p negative), the method stops without concluding. 

Second step. For each family, compute the indicial equation ( p. 36 ) and re- 
quire the necessary conditions : 

• C2. Every zero of detP (a Fuchs index) is integer. 

• C3. Every zero i of det P has a multiplicity equal to the dimension of the 
kernel of det P(i) 

V index i : (multiplicity of i) = dim Ker P(i). (5.113) 

Remark. There is no such condition as CI on page [ti], i. e. the indicial 
polynomial may have a degree smaller than N . If the indicial equation has 
degree N, the conditions C2 and C3 (N distinct integers in the one-dimensional 
case) are slightly stronger than the conditions in Bureau (JV— 1 distinct integers). 

The next step is easily computerizable 27, 4?J if one represents u^, . . ., 
as Laurent series bounded from below : u^^with powers in the range (p : +oo), 
with powers in the range (p + p : +oo), where p denotes the smallest Fuchs 
index, an integer lower than or equal to —1, ... 

Order n = is identical to the method of pole-like expansions and the 
Laurent series for u(°) 

u(o) =E u r^' +p ' ( 5 - m ) 

3=0 

represents a particular solution containing a number of arbitrary coefficients 
equal to one (index —1) plus the number of positive Fuchs indices, counting 
their multiplicity. 

Order n — 1 is identical to the "equation auxiliaire" of Darboux 

E« eE'(i,u("')u (1) -0, (5.115) 

and the Laurent series for 

uW =^ufV +p , (5.116) 

3=P 

represents a particular solution containing a number of arbitrary coefficients 
equal to the number of Fuchs indices, counting their multiplicity. If detP(«) 
has degree N, it represents the general solution of (5.115). 
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Consequently the sum u^ ) -feuW is already, in the neighborhood of (x, s) = 
(0, 0), a local representation of the greatest particular solution of ( 2.30| ) available 
in this method (the general solution if detP(i) has degree N), and this is a 
Laurent series with a strictly larger extension (p + p : +00) than that for the 
unperturbed expansion (p : +00). 

At each order n > 2, the sin gulari ty order of the particular solution of the 
linear inhomogeneous equation ( 5.110 ) E^ n > = is dictated by the contribution 
of the previously computed coefficients : it is increased by p at each order 
n and is equal to np + p 



+00 



Vn > : = uJ n V +p - (5.117) 



j=np 



(0,0) 



Third step. Solve the recurrence relation for U"- for all values of (n, j) 7^ 



Vn > Vj > np, (n,j) ? (0, 0) : E< n > = P(j)uf + Q™(x, {u^}) = 0. 



The generated necessary stability conditions are 
• C4. 



(5.118) 



Vn > V index i, (n, i) 7^ (0, 0) : Qj orthogonal to Ker adj P(i). 

(5.119) 

These orthogonality conditions must be satisfied whatever be the previ- 
ously introduced arbitrary coefficients. For a single equation, the condition 
C4 is simply 3 = 0. 

Qj"^ depends on all u^" ' with n' < n,j' — n'p < j — np, (n',j') 7^ (n,j), 
and this is the only ordering to be respected during the resolution. The costless 
ordering on (n, j) is the one which generates stability conditions the sooner, and 
it depends on the structure of indices of the DE under study. 

In order to avoid introducing more arbitrary coefficients than N, the precise 
rule is : 

• if n = or (n = 1 and i < 0), assign arbitrary values to mult(i) compo- 
nents of u^" 3 defining a basis of Ker P(i), 

• if (n = 1 and i > 0) or n > 2, assign the value to mult(i) components of 

defining a basis of Ker P(i). 

The resulting double expansion (Taylor in e, Laurent in x at each order in 
s) can be rewritten as a Laurent series in x extending to both infinities 



Vn > : u« = £ uJ n V +p , E W = 5] Ef V +q , (5-120) 

j=rip j=np 
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n=0 



+ 00 

fnl 

J = Tip 



J2 %X J+p . (5.121) 



j=-oo 



1. The Fuchsian perturbative method (as well as the nonFuchsian one which 



will be seen section 5^ ) is useful if and only if the zeroth order n — fails 
to describe the general solution. This may happen for two reasons. The 
most common one is a negative Fuchs index in addition to —1 counted 
once, the second, less common one is a multiplicity higher than one for 
some family, as in the example of section 5.7.3| . 



2. We do not know of an upper bound for n, but there exists a lower bound. 



Indeed, in the linear inhomogeneous ODE (5.11C), logarithms can arise 
only when some precise powers of x, only depending on the homogeneous 
part, are present in the rhs Laurent series R^ n \ The lower bound n results 
from the condition that the lowest Fuchs index and the highest one, once 
forced to interfere by the nonlinear terms, start to contribute to such 
dang erous powers. An example of such a condition is given in section 
5jJ. 



Even if all Fuchs indices are positive (except —1 counted once), the lower 
bound on n may be greater than 0, as in the example of section [5.7.3 . 



3. Remark on index — 1 . In the case of a single equation, since indices must be 
distinct integers, the condition Qp 1 ' = at the smallest Fuchs index i = p 
is identically satisfied. Nevertheless, the frequently encountered statement 
"resonance —1 is always compatible" is erroneous, and numerous nonzero 
stability conditions Q^l\ = can be found in the examples of f3lfl . Indeed, 
even at first perturbation order, the stability condition at index —1 may 
not be satisfied : just like Fuchs index i — p provides an identically satisfied 
stability condition, Painleve "resonance" —1 has the same property if and 
only if p = — 1, i.e. if — 1 is the smallest integer index. 

If p is different from — 1, Painleve resonance —1 seems to be satisfied, but 
it is only because a Laurent series ranging from power p to +00 cannot 
represent the general solution, thus preventing the building of Painleve 
stability condition Q^l = 0. 

5.7.1 Fuchs indices, Painleve "resonances" or Kowalevski 
exponents? 

Given a nonlinear algebraic DE of order N, one can define three sets of at 
most N numbers associated to it : 



1. the Fuchs indices of the auxiliary equation of Darboux (section 5.7) 
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2. the "resonances" of the nonlinear equation (section 5.4) 



3. the Kowalevski exponents, only defined if the nonlinear equation in invari- 
ant under a scaling transformation (x, u, E) — > (kX, k p u, fc q E). 

We have already seen the identity of the first two sets, defined for each family 
of movable singularities. Let us prove that the third notion is not distinct. The 
third set is defined as follows (for an introduction, see (l(J). The invariance 
implies the particular solution ("scaling solution") u' ' = const (x — xo) p , which 
is identical to a family of movable singularities. The Kowalevski exponents p 
are defined as the characteristic exponents of the linearized system near this 
solution, which proves the identity of the three notions. 

Said differently, all these numbers are Fuchs indices, and this link to the the- 
ory of linear DEs proves the uselessness of the notions of Kowalevski exponents 
and Painleve resonances. 

5.7.2 Understanding negative Fuchs indices 

The ODE with a meromorphic general solution |H| 

E = u xx + 3uu x + u 3 = 0, u = 1 -, a and b arbitrary, (5.122) 

x — a x — b 

has two families, 
(Fl) p = -1, 4° } = 1, indices (-1, 1), 
(F2) p = -1, 4 0) = 2, indices (-2,-1), 

and this provides a clear comprehension of negative Fuchs indices, since the 
index —2 must coexist with the meromorphy. Indeed, the representation of the 
general solution ( 5.122j ) as a Laurent series of x — xq is the sum of two copies 



of an expansion of l/(x — c), and there exist two expansions of l/(x — c) 



-l 



(x-c)- 1 = (c-xo)- 1 - J (x-x ) J , |c-x | < |x-x | (5.123) 

j=— OO 

= ^-( C -x ) 1 -J(x-x )^, |x-x | < |c-x |. (5.124) 

3=0 

The family (Fl) corresponds to the sum (first expansion with c = a — Xq) plus 
(second expansion c = b), while the family (F2) corresponds to the sum (first 
expansion with c = a) plus (second expansion with c = b). This can be checked 
by a direct application of the algorithm, which for family (F2) gives |3lj 
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+£ 3 { ^ x -7 + ^l x -6 + AlB 2 x - S 1^-4) + 0(£ 4 } (5 125) 

+ ei? lX - 2 + £ A lX - 3 + (^e 2 A 1 B 1 - ^B 3 ^ + 0(xX^6) 
2x - £Bl (5.127) 



X*-EB lX +$(-eA 1 + e*Bf)' 

where A\ and B\ are the arbitrary coefficients at order one. The simple pole 
X = with residue 2 has been "unfolded" by the perturbation into two simple 



poles with residue 1, at the two arbitrary locations | eB\ ± ^/2eA\ — e 1 B\ 
both close to 0. 

For other examples, see f9(| and conference proceedings referenced in plf . 

5.7.3 The simplest constructive example 

The equation 

u" + Auu' + 2u 3 = (5.128) 
is the simplest constructive example, because 

1. there exists a movable logarithm, as shown by the a— method (BSMF §13, 
P 221), 

2. the method of pole-like expansions fails to find it, 



3. the assumption of a "descending" Laurent series (5.34) fails to find it, 

4. the Fuchsian perturbative method finds it after a very short computation, 
as we now show. 

There exists a single family 

p = -1, = 4°Vo 0) - l) 2 = 0, indices (-1,0), (5.129) 

with the puzzling fact that Uq should be at the same time equal to 1 accord- 
ing to the equation Eq ^ — 0, and arbitrary according to the index 0. The 
application of the method provides 

i/ ' = x 1 (the series terminates) (5.130) 

E'{x,u {0) ) = dl+A X - 1 d x + 2 X - 2 (5.131) 

= u^x" 1 , arbitrary, (5.132) 

E^ = E'(x,uW)uW + 6n<%« 2 + 4u^u^' 

= X~ 2 (X 2 " (2) )" + 2« (0) V 3 = (5.133) 

u (2) - -24 1)2 ( X - 1 Log X - X - 1 ). (5.134) 
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The movable logarithmic branch point is therefore detected in a systematic way 
at order n = 2 and index i = 0. 

The necessity to perform a perturbation arises from the multiple root of the 
equation for , responsible for the insufficient number of arbitrary parameters 
in the zeroth order series . 



5.7.4 The two examples 



Example 1. "Complete (PI)" eq. ( 5.2C ). The met hod i s useless. 

Example 2 (Chazy complete equation of class III ( |5.24|) ). 

For the second family (in case 0I3 = 0), the method is useless. 

For the first family, since all indices are negative, one must start the pertur- 
bation process n > 1. 

To obtain the stability conditions up to a given order n > 1, we only need 
to compute the first 3n coefficients of each element : 

-3r+3ra-l 

«W= J2 «5'V~V = 0,... > n > (5.135) 

j=-3r 

,.(0) (1) (n-1) (n) ffil<un 

u 0:3n-l> "-3:3n-4J • ■ ■ > u -3n+3:2' u -3n:-2) ^d.lOOj 

where j\ : j'2 denotes a range of j values. The most efficient way to perform the 
double loop on (n,j) is to perform the outside loop in the variable k = j — np, 
with p = —3, and the precise double loop is : for k = to fcmax do for 
n = (if k = then 1 else 0) to rimax do solve the linear algebraic equation 
( |5.118D for uf\j = k + np. 

Let us compute all the stability conditions at first and second order. The 
computer printout reads (full details are given in Ref. ]3l[) : 

k = : Q a ] = 0, u W 3 arbitrary (5.137) 

k = 1 : Q W 2 ee -6(5ci + A2d z )u^\ = 0, arbitrary (5.138) 

k = 2 : g a l ee -12(ci + 9d 3 )M_2 + 18 (ci - 8d 3 )tt^' 
6 



(2c - 3cj - 117cid 3 - 594d§ + 18^ + 108d 3 )u a 3 = 05-139) 

o 

u^i arbitrary, (5.140) 
k = 3 : ee _-^ lU W 2 = 0, (5.141) 





k = 4 : ee -(8do + 57d[)u°:l ~ Udiu^u^l = (5.142) 
1 m 2 24 



fc = 5: ee - — {18d' + 99d'()u°:l - —dm^l 

o 

+ — (16do + 72d' 1 )wL 1 ]uL 1 2 = 0. (5.143) 
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Five conditions are obtained, three at order one, equivalent to (I3 = c\ = 
co = 0, and two at order two, equivalent to d\ = do = [in order to simplify 
expressions, we have put the first order conditions in the above expressions for 
k = j + in > 3], after seventeen values of (n,j). These conditions were given 
without any det ail by Chazy p^L T hey restrict the complete ODE (5.24) to the 
simplified ODE ( |l.l7[ ), modulo (||]). 



Chazy proved the general solution of ( 1.17 ) to be single valued inside or 
outside a circle whose centre and radius depend on the choice of the three 
arbitrary constants; it is holomorphic in this domain, and the only singularity 
is a movable natural boundary ("coupure essentielle" ) defined by this circle. He 
also gave a parametric representation of the general solution u(x) in terms of 
two solutions of the (linear) hypergeometric equation, but single valuedness is 
not at all apparent on this representation. 

The direct explicit solution of Bureau 0, [ll| is given section 5.5.1 . 



5.7.5 An example needing order seven to conclude 

The following equation, isolated by Bureau ([[[B) p. 79), 

u"" + 3uu" - 4u' 2 = (5.144) 

possesses the two families 

p = -2, u ( 0) = -60, ind. (-3, -2,-1, 20), K = u"" + 3uu" - 4^8,145) 

(5.146) 



p = —3, Uq ^ arbitrary, indices (—1, 0), K = iuu" — \v! A . 
The second family has a Laurent series (p : +00) which happens to terminate 



m 



i' ' = c(x — xq) 3 — 60(a; — xa) 2 , (c, xo) arbitrary. 



(5.147) 



The Fuchsian perturbative method is useless, for the two arbitrary coefficients 
corresponding to the two Fuchs indices are already present at zeroth order. 

The first family provides, at zeroth order, only a two-parameter expansion 
and, when one checks the existence of the perturbed solution 



+00 

E« 

71=0 



+ OO 



3=0 



(5.148) 



one finds that coefficients u^ \ vP^, vS^, u_\ can be chosen arbitrarily, and, at 
order n — 7, one finds two violations pi| 



n( 7 ) = f# (°> w ( l > - n 

V-l — "20 -3 — u ' ^20 



(7) = u (orji)\,(i) _ 



i 20 U -3 U -2 — 0, 

implying the existence of a movable logarithmic branch point. 



(5.149) 
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Remark. The value n = 7 is the root of the linear equation n(i m i n — p) + 
(imax — p) = — 1, with p = -2, i m in = -5, i ma x = 18, linking the pole order p in 
the Fuchsian case c = 0, the smallest and the greatest Fuchs indices. It expresses 
the condition for the first occurence of a power x , leading by integration to 
a logarithm, in the r.h.s. R^ of the linear inhomogeneous equation (5.110), 
r.h.s. created by the nonlinear terms 3uu" — 4u' 2 . 



5.7.6 Closed-form solutions of the Bianchi IX model 

In this example, the no-log conditions are used in a constructive way, in 
order to isolate all possible single valued solutions. 

The Bianchi IX cosmological model fFq] is a system of three second order 
ODEs 

(Log A)" = A 2 -{B- C) 2 and cyclically, ' = d/dr, (5.150) 
or equivalently 

/T \ll 2,2 2 2 / 2 /1 / 2 

(LogWi) = 0J 2 + ^3 — ^2^3/^1 J A = GJ2U13/0JI, W\ 

One of the families g, |7| 

^ = X _1 +a 2X + 0(x 3 ), X = t 
B = X" 1 +&2X + 0(x 3 ), 
C = X" 1 +C2X + 0(x 3 ), 

has the Fuchs indices (— 1, —1, —1, 2, 2, 2). The Fuchsian perturbative method 

N 2+N-n 

A = x~ 1 ' S ^2,£ n a ^X^ X = r ~ T 2, and cyclically, (5.153) 

n— j ——n 

then gives a failure of condition C4 at (n, i) = (3, —1) and (5, —1) [^9|, and the 
satisfaction of these no-log conditions generates the three solutions : 

(6<°) = 4 0) and b a l = c^\) or cyclically (5.154) 
4°) = 4 0) = 4° } - 0, (5.155) 
W =6 W =C W (5.156) 

These are constraints which reduce the number of arbitrary coefficients to, re- 
spectively, four, three and four, thus defining particular solutions which may 
have no movable critical points. The question is : do they define additional 
solutions to what is known? 

The only three closed-form solutions which are known are single valued, they 
are defined as the general solution of the following three subsystems. 



= BC and cyclically. 

(5.151) 



72, 

(5.152) 
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1. The 4-dim axisymmetric case B = C 10S], whose general solution (2.60) 
is trigonometric. 



2. The 3-dim Darboux-Halphen system |44| , pl|] 

lu[ = lu2^3 — tjJ\L02 — ^i^3 7 and cyclically. 



(5.157) 



3. The 3-dim Euler system (1750) 
around its center of mass 



describing the motion of a rigid body 



uj[ = W2UJ3, and cyclically, 



(5.158) 



whose general solution is elliptic (8), see (2.31) and ( 2.45| ). 

The first constraint ( [5.154 ) implies the equality of two of the components 
( A, B , C) at every order and thus represents the four-parameter solution of Taub 



The second constraint ( 5.155) represents the three-parameter solution of the 

Darboux-Halphen system ( 5.157 ). 

The third and last constraint ( 5.156| ) represents an extrapolation to four 
parameters of the three-parameter solution ot the Euler system described by 
= 0. This would-be four-parameter, global, closed form, single 



,(o) + 6 (o) 



valued exact solution has not yet been found. 



5.8 The nonFuchsian perturbative method 

Whenever the family under study has a number of Fuchs indices smaller 
than the differential order N, the Fuchsian perturbation method fails to build 
a representation of the general solution, thus possibly missing some stability 
conditions. Examples are ( 5.146) and the second family of ( |5.24| ) in the case 



0. The missing solutions of the auxiliary equation ( 5.115| ) are then nonFuchsian 
solutions, see section 5.2.2) . 



There is no difficulty to algorithmically compute the nonFuchsian expansions 
( |5.11 ), but these are of no immediate help, due to their generic divergence. 



There is one situation where some stability conditions can be generated 
algorithmically (indeed, we are not interested in computations adapted to a given 
equation, only in computerizable methods). It occurs when the two following 
conditions are met | p5[ . 

1. There exists a particular solution u = u*- - 1 which is known globally, mero- 
morphic and has at least one movable pole at a finite distance denoted 
x . 

2. The only singular points of the linearized equation E^ 1 -* = are x — Xo, 
nonFuchsian, and x = 00, Fuchsian. 
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Then, the property that a fundamental set of solutions be locally single 
valued near x — x — x o = is equivalent to the same property near x = 00 • 
This is the global nature of u*- ' which allows the study of the point \ = 00, 
easy to perform with the Fuchsian perturbation method. 

An important technical bonus is the lowering of the differential order N of 
equation E^ — by the number M of arbitrary parameters c which appear in 
u(°). Indeed, again since u*- - 1 is closed form, its partial derivatives 9 c u(°) are 
closed form and are particular solutions of E^ 1 ' =0, which allows this lowering 
of the order. 

At each higher perturbation order n > 2, one similarly builds particular 
solutions u.(") as expansions near \ — 00 an d one requires the same properties. 



5.8.1 An explanatory example : Chazy's class III (N 

3,M = 2) 



The simplified equation (1.17), which possesses the PP p2j and therefore for 
which no is multivalued, is quite useful just to un derst and the method. This 
equation admits the global two-parameter solution ( 5.57 ) = c\~ 2 — &X 1 - 
The linearized equation 

= E'{xM 0) )u {1) = [dl - 2u^d 2 x +6u^d x - 2u^]u^ = (5.159) 

possesses the two single valued global solutions d X(s u^ and d c u^ \ i.e. vS 1 ' — 
X~ 3 ,X~ 2 j an d it has only two singular points \ — (Fuchsian) and x = 00 
(nonFuchsian with Thome rank two). The lowering by M = 2 units of the order 
of the linearized equation results from the change of function 



a 



« = X ~ 3 v ■ = X 3 [d, + - 2c X - 2 ]v" = 0, (5.160) 



and the study of the Fuchsian point x = 00 yields an integer Fuchs index, which 
proves the global single valuedness of the general solution u^K 
Remarks. 



• The local study of x — provides a formal expansion (5.11) which happens 
to terminate, a nongeneric situation, thus providing the fundamental set 
of global solutions at perturbation order n = 1 

V X Vc : u« = X ~ 2 , X" 3 , {e~ 2c/x - 1 + 2c X ' 1 )x- 2 1 '(2c 2 X-5.161) 

This proves the existence of an essential singularity at x — ( re f- [ |66[ 
chap. XVII). 

• Going on with the formalism of Painleve's lemma at higher orders consti- 
tutes the rigorous mathematical framework of the local representation of 
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the general solution obtained by Joshi and Kruskal |7Q] 

" - -**r l +<** (' + • - j + m - 7 t4 + 0(E,) ) ■ « - ? 2 * 

(5.162) 

This representation reduces to the one given by Chazy (Taylor series in 
1/x) if one starts from the Fuchsian family u ~ — 6x _1 - 

5.8.2 The fourth order equation of Bureau (N = 4, M = 2) 



In section 5.7.5 , the fourth order equation ( 5.144 ) has been proven to be 
unstable after a computation practically untractable without a computer. Let 
us now prove this r esult w ithout computation at all Q. For the global two- 
parameter solution (5.147), the linearized equation 

= E'(x, uW)uW = [fig + 3u(°>fi£ - 8u^d x + 3u®>] u W = (5.163) 

has only two singular points x = (nonFuchsian) and x = oo (Fuchsian), it 
admits the two global single valued solutions d Xo u^ and d c u^\ i.e. uW — 
X~ 4 , X -3 - The lowering by M = 2 units of the order of the linearized equation 
( [5.163 ) is obtained with 



= x 4 v : 



[d 2 x - l&x^d x + 3c X - 3 - 60 X - 2 ]«" = 0, 



(5.164) 



and the local study of x = °o is unnecessary, since one recognizes the confluent 
hypergeometric equation. The two other solutions in global form are 







= x- 3 oF 1 (24;-3c/x) = X 17/2 J23(V^/x), (5-165) 
= X 17/2 N 23 (^/12^), (5.166) 



where the hypergeometric fonction o-Fi(24; — 3c/%) is single valued and possesses 
an isolated essential singularity at x = Q, while the fonction N23 of Neumann is 
multivalued because of a Log x t erm. 

Re mark . The local study of ( 5.163 ) near x — provides the formal expan- 
sions ( 5.11 ) for the two nonFuchsian solutions 



X - 0, c ? : n« - eW-iVx//^! + (^)) ; 



(5.167) 



detecting the presence in ( |5.163 ) of an essential singularity at x = 0, but the 
generically null radius of convergence of the formal series forbids to conclude 
to the multivaluedness of A nonobvious result is the ex istence , as seen 

above, of a linear combination of the two formal expansions ( [5.167 ) which is 
single valued. 
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5.8.3 An example in cosmology : Bianchi IX (N — 6, M 

4) 



The Bianchi IX cosmological model in vacuum ( 5.15C) d oes not possess the 



Taub Jl08| found the general 



PP [gjj |7J. Let us prove it rapidly 
solution of the axisymmetric case of two equal components, a meromorphic 
expression ( 2.60| ) depending on the four arbitrary parameters (k%, fca, t\, T2). 
The linearized system generated by the perturbation 



A = A {0) (1 + eA (1) + 0{e 2 )) and cyclically 



(5.168) 



has the differential order N = 6, which is then lowered by M — 4 units by the 
change of function dictated by the symmetry of the system : PW = B^ + 
CW.MM = pW-C« 

A^"~2A^ 2 A^ = 0, (5.169) 
P (D" _ 2A^B^P^ = 4(A^B^ - A^')A^\ (5.170) 
M W" + 2(#5(°) - 2P( 0)2 )M (1) = 0. (5.171) 

Indeed, the four single valued global solutions 

(AW ) pW)=S c (Log^ ),Log(p(°) +C<°>)), c = A 1 ,* a ,ri,75, (5.172) 

are those of the equations ( 5.169| )-( 5.170| ), 



Aff 1 ) -0,(A( 1 ),P( 1 ) + 2A( 1 )) 



((t - ri) cothfci(r - n) - 1/fci, 0), 
(0, (t - r 2 ) cothfc 2 (r - t 2 ) - l/fc 2 ), 
(cothfci(r — ri),0), 
(0, cothfc 2 (r - t 2 )), 

(5.173) 

and there only remains to study the equation (5.171). It has a countable infinity 
of singular points : r — r 2 = im7r/fc 2 , m£Z (nonFuchsian, of Thome rank two), 
accumulating at r = 00. This uneasy situation can be overcome by taking the 
limit ki = fc 2 = 0; it is indeed sufficient to exhibit a movable logarithm in this 
limit, for it will persist for (ki, fc 2 ) ^ (0, 0). In this limit 



h = k 2 = : 



At 2 



4ft -1) 



M« = 0, t = — %.174) 

n - t 2 



the only singular points are t = (nonFuchsian) and £ = 00 (Fuchsian), the 
optimal situation. The Fuchs indices being —2 and 1, the computation of three 
terms is sufficient to exhibit a logarithm, and this proves the absence of the 
Painleve property for the Bianchi IX model in vacuum. 

Remarks 



77 



1. The two solutions are globally known |?{| : 

h = k 2 = : M« = e- 2 ' l r\ e^H- 1 t ' z- 4 e 4z dz(5.175) 

which shows the presence of a logarithmic branch point at t = 0, or at 
t = oo as well. 

2. The two formal non-Fuchsian solutions are 

+oo 

t - t 2 -» : M« - e Q /(--^) ^ A fc (r - r 2 )' £ + s , A ^ 0, (5.176) 

fe=0 

with 

a = ±2£;f 1 sinh£;i(T2 - n), s = 1 =F 2coshfci(r 2 - n). (5.177) 

The two generically irrational values for the Thome exponents s allow to 
conclude only if the divergent series Afe(r — r 2 ) k can be summed. 



5.9 Miscellaneous perturbations 

The differential complexity of the a— method explains why it usually suc- 
ceeds in case of failure of all the other methods, which only have an algebraic 
complexity. Consider the ODEs, none of which admits a power-law leading 
behaviour 

-2uu" + 3u' 2 + d 3 u 3 = 0, d 3 ^ 0, (5.178) 
u'" + uu" - 2u' 2 = 0, (5.179) 

(5.180) 
The first one is 

extracted from Chazy's class III (5.24) by the perturbation u = e~ 1 U 1 x = 
xq + eX, and it represents its second family, see section 5.4.1. The second and 



III! in II o n n 

u + Zuu — iu = U, 
and let us prove that each of them has movable logarithms. 



third ones were considered by Chazy [|23j, [24| who had to establish a special 
theorem, using divergent series, to exhibit the movable logarithms. Having 
degree one, none of th ese OD Es admits singular solutions. 

The first equation ( 5.1 T8j ) is classically processed by the a— method 



(n) 



E : 



n=0 



71=0 



■ n E^ n \ x = x +sX, (5.181) 



resulting in 

E (o) 

E (i) 



c(X-X a )- 2 , 



3u^' 2 = 
(Xq,c) arbitrary, 



(5.182) 
(5.183) 



= c(X - X )- 5 [-2((X - Xofu^Y' + c 2 d 3 /(X - X )\ = (#.184) 
= c 2 d 3 (X-X ) 3 [(X-X )Log(X-X Q )-{X-X )}/2. (5.185) 
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and proving the instability at perturbation order one. 



For equations (5.179) and (5.180), there exists no perturbation satisfying 
the assumptions of Theorem II page EE), there only exist singular perturbations, 
i. e. which discard the highest derivative. Since they however give the correct 
information, i t woul d be desirable to extend Theorem II in that direction. 

Equation ( [5.179] ) is handled by the singular perturbation 



+00 



+00 

E = e- 2 £ ™^ (r ' 



(5.186) 



which excludes u'" from the simplified equation 



E 



(0) 



ur> - 2u (0) ' 2 = 

X = x — xq, (xq, c) arbitrary, 
= c( X - 3 ( X V 1 ))"-6 X - 4 ) = 0, 

„(i) 



« (0) = c X -\ 



6x X (Logx- 1) 



(5.187) 
(5.188) 
(5.189) 
(5.190) 



This same perturbation (5.186) solves the case of the equation ( 5.18G| ) 

E^ ee 2u^u m " - 3u^' 2 = (5.191) 

(5.192) 

= (5.193) 
(5.194) 



,.(°) = 



cx 



-60 X ~ 2 (Logx-l). 



5.10 The perturbation of the continuum limit of 
a discrete equation 

Discrete equations can be considered as functional equations linking the val- 
ues taken by some field variable u at a finite number N + 1 of points, either 
arithmetically consecutive : x+kh, fc— fco = 0, 1, . . . , N, or geometrically consec- 
utive : xq k ,k — k = 0, 1, . . . , N, where h or q is the lattice stepsize, assumed to 
lay in some neighborhood of, respectively, or 1, and fco is just some convenient 
origin. 

Definition [[35|. A discrete equation is said to possess the discrete Painleve 
property if and only if there exists a neighborhood of h = at every point of 
which the general solution x — > u(x,h) has no movable critical singularities. 

Consider an arbitrary discrete equation ( |5.195 ), 

MxMh : E(x,h,{u(x + kh), k - k = 0, . . . , N}) = (5.195) 

algebraic in the values of the field variable, with coefficients analytic in x, the 
stepsize and some parameters a. Let (x, h, u, a) — * (X,H,U,A,s) be an ar- 
bitrary perturbation admissible by the suitable extension of the theorem of 
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Poincare to discrete systems. Two such perturbations are well known, the au- 
tonomous limit 

x = xq + sX, h = eH, u=U, a= analytic (A, e), (5.196) 
and the continuum limit 

x unchanged, h = s, u = U, a= analytic (A,e). (5.197) 

The latter can be extended into a perturbation of the continuum limit p5| 

+00 

ir unchanged, h = e, u = e n u^ n \ a — analytic (A, s), (5.198) 



entirely analogous to the Fuchsian (section |5.7| ) or nonFuchsian (section 
perturbative method of the continuous case. 

It generates an infinite sequence of (continuous) differential equations = 
whose first one n = is the continuum limit. The next ones n > 1, which are 
linear inhomogeneous, have the same homogeneous part 

= indepen- 
dent of n, defined by the derivative of the equation of the continuum limit, while 
their inhomogeneous part ("right-hand side") comes at the same time from 
the nonlinearities and the discretization. 

Let us just handle the Euler scheme for the Bernoulli equation 

E = (u- u)/h + u 2 = (5.199) 

(notation is u = u(x) 1 u = u(x + h)), i.e. the logistic ma p of Vc rhulst, a paradigm 
of chaotic be haviou r. Let us expand the terms of ( 5.199 ) according to the 
perturbation (5.198) up to an order in e sufficient to build the first equation 
E^> = beyond the continuum limit E^ = 

u = m (0) +M (1) e + 0(e 2 ) (5.200) 

u 2 = u^ 2 +2u^u^e + 0(e 2 ) (5.201) 

u = u + u'h+ (l/2)u"h 2 + 0(h 3 ) (5.202) 



= u (oy 



+ (u w + (1/2)«<°)")e + 0{e 2 ). (5.203) 



h 

The equations of orders n = et n = 1 are written as 

= M (°)'+ u (°) 2 = (5.204) 

£(D = £(°)' u (i) + (l/2) u (°)" = 0, E^' =d x + 2«W. (5.205) 
Their general solution is 

,X = x ~ x 0i x o arbitrary (5.206) 



u (°) = X 



11 



(i) = J 1 



u_{x ~X Log"0, ip — x — xq, u_{ arbitrary, (5.207) 



and the movable logarithm proves the instability as soon as order n — 1, at the 
Fuchs index i — — 1. 
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5.11 The diophantine conditions 

In fulfilling the systematic programme of Painleve, one encounters the fol- 
lowing kind of diophantine equations 



£- = i, (5.208) 



fe=l 

with n and p given integers, whose unknowns (n^) are Fuchs indices, which 
must therefore be either integer or infinite. They admit a finite set of solutions, 
which allows all cases to be further examined. Details can be found in Bureau 
1964, M.I, MIL 

Such a diophantine condition always arises when there exists more than one 
family, as the constraint that, simultaneously, all Fuchs indices of all families 
be integer. Let us just give one example pi]| . The Hamiltonian Henon-Hciles 
system j32j in two coupled variables (91,92) 

H = (l/2)(5? ia + q\ x + C 1 q 2 +C 2 q 2 ) + aq 1 q 2 -(l/3)Pq 3 1 =E, (5.209) 
qi,xx + c x qx - (3q\ + aqj = (5.210) 
q2,xx + c 2 q 2 + 2aqiq 2 = (5.211) 

defines by elimination the fourth order ODE in v = q\ J5(J 

Vxxxx + (8a - 2/3)vv xx - 2(q + p)v 2 x - (20/3)a^ 3 

+(ci + 4:C 2 )v liXX + (6aci - A[3c 2 )v 2 + Adc 2 v + 4aE = 0, (5.212) 

with (a, /?, ci, C2, E) constants. Let us restrict here to c\ — c 2 = 0. 
There exist two families 

3 

p = —l,Vo = —, indices (— 1, 10, n, r 2 ), (5.213) 
a 

6 

p = -1, w = -— indices (-1, 5, Si, s 2 ), (5.214) 
P, 

in which and Si satisfy the equations 

r 2 - 5r + 12 + 67 = 0, s 2 - 10s + 24 + 487" 1 = 0, 7 = 0/a. (5.215) 

The diophantine equations to be solved are 

(ri - r 2 ) 2 - (2fc - l) 2 , (si - s 2 ) 2 = (2Z) 2 , (r» + l)(r, - 10)(s 4 + l)( Si - 5) ^ 0, 

(5.216) 



r 



with fc and i two strictly positive integers. Making use of (ri — r 2 ) = (2 
5) 2 , (si — s 2 ) 2 = 4(s — 5) 2 , the elimination of 7 between ( 5.215 ) yields 

48 , 1152 

^=t^p> l =1+ 23 +i 2k-ir (5 - 217) 
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and this provides sharp bounds for I : 1 < I 2 < 49. One thus obtains the four 
solutions for j3/a, (fc, I), (n, 7-2), (si, S2) 



- 1 
-2 
-6 
-16 



(1,7), (2,3), (-2,12) (SK), 
(3,5), (0,5), (0,10), 
(6,3), (-3,8), (2,8) (KdV5), 
(10,2), (-7,12), (3,7) (KK). 



(5.218) 
(5.219) 
(5.220) 
(5.221) 



Three of them restrict the ODE to the stationary reduction of w ell-known 
soliton equations, thus proving the PP : Sawada-Kotera (SK |103|), higher- 
order Korteweg-de Vries (KdV5, f(|) and Kaup-Kupershmidt (KK, @ |l|) 
equations. 

The case j3 = — 2a is similar to that of the ODE in section 5.7.3 : «o is a 
double root of its algebraic equation and is not arbitrary although is an index. 
The results of the Fuchsian pcrturbative method are also similar; listed by 
increasing cost (number of needed values of (n,i)), the first stability conditions 
q| ti) = are 



(i) = 



Qo 

o {0) 



0, cost = 2 
-40a4 1)2 = 



0, cost 



-30a 3 4 0) = 0, cost = 10 



T20au 



(0) (l) 



0, cost = 12. 



(5.222) 
(5.223) 
(5.224) 
(5.225) 



To detect the instability, the method of pole-like expansions is here sufficient 
but the Fuchsian perturbative method is much cheaper. 
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Chapter 6 

Construction of necessary 
conditions. The Painleve 
test 



This chapter makes the synthesis of all the methods of chapter |^ in order to 
define a usable end product which makes obsolete the meromorphy test of section 
|2.3| . This end product is widely known as the Painleve test. Before detailing 
the steps of this algorithm in section 6.6, for ODEs as well as for PDEs 



some 

implementation of physicists' 



prerequisite technical developments are needed 
desiderata (section 6.1), technicalities to simplify the computations (section 6.2) 
and the qui te important feature of the invariant Painleve analysis (sections 6.3, 
|Jand|6l|. 



6.1 Physical considerations 



Some DEs encountered in physics are unstable, although integrable or par- 
tially integrable in some obvious physical sense. It is then extremely important 
not to discard them; this is achieved by relaxing some of the mathematical 
requirements. 

Firstly, nonpolynomial DEs can be made polynomial by transformations on 
u like in section 



3.3.3 , neces sarily outside the groups of invariance of the PP 
defined in sections 3.3.1 and |3.3.2 . 



Example 1 (sine-Gordon). 



(sine-Gordon) 



sin it, e 



2(w xt - v x v t ) - v 3 + ' 



(6.1) 
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Example 2 (Bcnjamin-Ono). The nonlocal, nonpolynomial PDE 

1 f +oc v(x' t) 
ut + uu x + U(u xx ) = 0, H(u) = -pp / -i-L^da/, (6.2) 

7T J_ OQ X'-X 

in which H is the Hilbcrt transform, and pp the Cauchy principal value distri- 



bution, is equivalent fcOL 101] to the local and polynomial system 



U t + UU X + U X y, U XX + Uyy = (6.3) 

in one additional independent variable y. 

Secondly, unstable polynomial DEs may be made stable and polynomial by 
transformations like fl3.9|). 

Example 3 (parity invariance). The Ermakov-Pinney ODE jl9], ^ 

u xx - a 2 u + /3 2 u~ 3 = 0. (6.4) 

is unstable (algebraic branch point p = 1/2) and invariant by parity on u : 
the transformation u — > u 2 or u~ 2 preserves its polynomial form and makes it 
stable. 



6.2 Technicalities 

A careful choice of the dependent variables can save many computations. 

Example 1 (dynamical systems). These systems of first order ODEs some- 
times possess an equivalent scalar ODE. This is the case of the Lorenz model 
(O), equivalent to Jl07 



xx"'-x'x"+x 3 x'+(b+(7+l)xx"+{a+l)(bxx'-x' 2 )+ax i +b{l-r)crx 2 = (6.5) 



and of the Henon-Heiles Hamiltonian system in (91,172) ( 5.209Q which implies 



the fourth order ODE in qi only (5.212). This offers two advantages. The first 



one is to reduce the matricial recurrence relation to a scalar one. The second one 
is much more interesting : the scalar ODE has a number of families lower than 
or equal to that of the DS, which saves a lot of useless cases to consider; thus, 
in the HH system, the leading powers for q 2 ) are (—2, —2), (—2, —1), (—2, 0), 
while the equivalent fourth order ODE for qi has only one leading power —2. 

Choosing an integrated dependent variable for the computations saves a lot. 
The principle is that, if a DE for u is to be stable, this allows the presence of 
one movable logarithm in its primitive v = J u dx. If changing u to v x allows 
the DE to be integrated once or more, expressions are shortened. 

6.3 Equivalence of three fundamental ODEs 

Let <S* be a given analytic function of a complex variable x, and let us consider 
the three differential equations in ip, \, ip : 
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^fxxx 


3 / fxx\ 2 


= s 




2 I <p x J 




-x-\ 


~2u> x — 2uj 2 


= s 




n ^XX 


= s 







(6.6) 
(6.7) 
(6.8) 

The first one is the Schwarz equation; if read backwards, it defines 5* as the 
Schwarzian {ip; x} of (p. The second one is the Riccati equation in its normalized 
form (equations for u> or \ are equivalent and both of Riccati type). The third 
one is the second order linear Sturm-Liouville ODE in its normalized form. 

Each of these three ODEs possesses a fundamental uniqueness property. 
As shown by S. Lie, the Schwarzian is the unique elementary homographic 
differential invariant of a function ip, i.e. the unique elementary function of the 
derivatives Dip of ip, excluding ip itself, invariant under the 6-parameter group 
Ti. (or Mobius group, or PSL(2,C)) of homographic transformations : 

Tt : ip — > — — — — -, (a, b, c, d) arbitrary complex constants, ad — be = 1. (6.9) 

cip + d 

Among nonlinear first order ODEs in the class : 

u' = R(u,x) (6.10) 

where R is rational in u and analytic in x, the Riccati equation is the unique 
one whose general integral has no movable critical points. As to equation ( |6.8| ), 
its uniqueness lies in its linear form. 

It is a classical result due to Painleve (1895, Legons p. 230, Oeuvres I) that 
the three ODEs fl6.6|), (6.7), ( |6.8[ ) are equivalent : it is sufficient to integrate 



one in order to integrate the two others. Consequently, any ODE reducible to 
one of these three ODEs can be considered as explicitly linearizable. The six 
ODEs obtained by elimination of S between any two of the three ODEs have 
the general solution : 

( \ Clfx tfxx I a -. -. x 
C!ip + C 2 2lf x 

V#) = (cip + cater* (6.12) 

<P(v) = x- ) r, cic 4 - c 2 c 3 = 1 (6.13) 

C 3 (w 2 - Wl) + C 4 (tJ 3 - UJl) 

1p{u) = Cllpl + 02^2, "01 = 7 ^ w — 7,^2=^1— —^6.14) 

{UI2 - Wi)(W3 - OJi) UOj,-L0 2 

VKV) = — ; — ; — , C1C4 - c 2 c 3 = 1 (6.15) 

uj(ip) = — (6.16) 

ClVl + C2^2 
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where the Ci 's are arbitrary constants, u>i and ipi particular solutions of (6.7) 
and (6.1 



Only two of these six solutions, namely x(v) an d , 4 ) { l p), eq. ( 5.1l| )-( |6.12 ), 
are expressed with a single function; therefore, among the three equivalent func- 
tions, <p is the most elementary one, and we are going to see that the two others, 
x(y) an d are t ne basic building blocks of the invariant Painleve analysis 

of both PDEs and ODEs. 

If the space of independent variables is multidimensional, for each additional 
independent variable t let us define a function C(x,t, . . .) by : 



(6.17) 



As seen from eq. ( 6.11 )— ( 3.12 ), the t dependence of the three equivalent functions 
is then characterized by the three equivalent linear PDEs (two homogeneous, 
one inhomogeneous) : 



ft 



■C<p x =0 
-l 



X 



(Cuj - -C x ) x = 



1 



4>t + Ctp x - -c x ip 



o. 



(6.18) 
(6.19) 

(6.20) 



The linearity of these PDEs, as well as the invariance of C under the change of 
function ip — > F(<p), F arbitrary, show that all independent variables but one 
give rise to linear equations. 

Systems (^)-(§]| and (|6~T8|)-(|6~20|) require the cross-derivative condition : 

Vx 1 ((fPxxx)t - (Vthxx) = 2((x _1 ) t )x - 2(( X ~ 1 ) x ) t (6.21) 
- 2i>- 1 ({i> t ) xx - {ip xx ) t ) =St + C xxx + 2C X S + CS X = 0. (6.22) 



6.4 Optimal choice of the expansion variable 

A PDE has movable singularities which are not isolated, on the contrary to 
an ODE, but which lay on a codimension one manifold 

ip(x,t,...)-ip =0, (6.23) 

in which ip is an arbitrary function of the independent variables and i^o an 
arbitrary movable constant. Even in the ODE case, the movable singularity can 
be defined as f(x) — (po, since the implicit functions theorem allows this to be 
inverted to x — xq — 0; this provides a gaude freedom to be used later on in 
chapter 0. 

The singular manifold and the expansion variable play two different roles, 
and there is no a priori reason to confuse them, so let us denote ip the function 
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which defines the movable singular manifold tp — tpo = 0, and x the expansion 
variable. The only requirement on \ is that it must vanish as tp — tpo and be a 
single valued function of tp — tpo and its derivatives. 
The Laurent series for u and E are defined as 



E 



3=0 



(6.24) 
(6.25) 



To illustrate our point, let us take as an example the Korteweg-de Vries 
equation 



E = 



+ 6uu x = 



(6.26) 



(this is one of the very rare locations where this equation can be taken as an 
example; indeed, usually, things work so nicely for KdV that it is hazardous to 
draw general conclusions from its s ingl e study). 

With the choice x = f ~ fo | 111 | , the coefficients (uj,Ej) are invariant 
under the two-parameter group of translations tp —f tp + b, b arbitrary complex 
constant, and therefore they only depend on the differential invariant ip x of this 
group and its derivatives : 



2—2,0 —lift ^ <Pxxi 

U=~2fxX +fyxxX +7. o 

6<p x 3 tp x 



(Pxx 



+ 0( X ). (6.27) 



F(ip), F arbitrary function, 



[The quantity C = is invariant under <p 

and therefore is uninteresting for the moment.] 

With the choice x — if ~ fo)/fx, always possible since the gradient of tp 
has at least one nonzero component, the invariance is extended to the four- 
parameter group of affinc transformations tp — > atp + b, (a, b) arbitrary complex 
constants, with accordingly a dependence on the differential invariant tp xx /ip x 
and its derivatives : 



tp x 6tp x 3 



ft 





1 




. fx _ 


*~6 


. fx _ 



+ 0( X ). (6.28) 



Let us extend this invariance to the six-parameter homographic group. 
Eliminating tpo between x and Xx for each of the two choices of x, one obtains 
the ODEs of order one for x 



Xx - Vx 

- Xx 



X 



^Pxx 
<Px 



o, x = 

= 0, x 



tp-tpo 
= f-fo 

<Px 



(6.29) 
(6.30) 
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whose coefficients only depend on the respective differential invariants. Now, 
one knows since S. Lie the differential invariant of the homographic group 



S = {p;x} = 



fxx 


1 


fxx 


. fx _ 


x~2 


. fx _ 



(6.31) 



and ipso facto the associated ODE of order one. Its general solution leads, 
by taking the homographic transform which vanishes as ip — ipo, to the good 
expansion variable 



f- fo 



fx 



l^if-fo) 



fx 



f-fo 



^XX 

2fx 



(6.32) 



Check : due to the homographic dependence of \ on ip, grad % is a polyno- 
mial of degree two in x with coefficients homographic invariants. Denoting t an 
arbitrary independent variable, possibly equal to x, one obtains 



1 



Xt 

Xx 

(Log ip) t 



-c + C xX 

z 

1 S o 

-Cx- 1 + \c x 



{CS + C xx ) x 2 



-C( Log ip) x + -C x 



(6.33) 
(6.34) 
(6.35) 



Again, eq. ( |6.33| )-( 3.34 ) arc not different from eq. (6/7), ( 6.19 ). 
The only price to pay for invariance is to privilege some coordinate x. 
For our KdV example, the final Laurent series, to be compared with the 
initial one (3.27), is remarkably simple : 



~ o C 25 1 _ , > 
-2x- 2 -g- T+ 0( X ). 



(6.36) 



The successive values of x an( i the corresponding subgroup items are gath- 
ered in the following table. 



Group 


Invariant I 


Riccati(x) 


Solution for \ 


<p + b 


fx 


Xx = I 


f-fo 


aip + b 


fxx/fx 


Xx = i - ix 


if - fa)/ fx 


(atp + b)/{c<p + d) 


{<p;x} 


Xx = l + (I/2) X 2 


<p — <po 



KruskaVs choice 

Kruskal |68) indicated the very simple choice x = x ~ • • •) °f expansion 
variable to make the practical computations as short as possible. This choice 
is equivalent in our formalism to a choice of gauge, namely S = 0, C' x = 0, 
and ip is then an arbitrary homographic function of x — fit, . . .) with constant 
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coefficients. The choice of Kruskal is really a choice of the expansion variable, 
not of the singular manifold, i.e. Xx = 1, not ip x = 1. 

Caution : this choice should only be used at the stage of building neces- 
sary conditions for the PP, and never at the stage of sufficiency because of the 
constraints put on (5, C) . 

6.5 Unified invariant Painleve analysis (ODEs, 
PDEs) 

This is a reference section containing all the items of that version of Painleve 
analysis which is common to ODEs and PDEs and which generates the simplest 
possible expressions, due to its built-in invariance. 

Consider a DE 

E(u,x) = (6.37) 

polynomial in u and its derivatives, analytic in x (E, u, x multidimensional), and 
the Laurent series for u and E around the movable singular manifold ip— tpo = : 

u = u_ p , 1 Log^ + ^u J x J+p J -peZ (6.38) 

3=0 

+ 00 

E = ^E jX j+q , -qeZ (6.39) 

3=0 

The coefficient u_ p i can be nonzero only if E does not explicitly depend on 
u. Let us denote x any independent variable such that (p x ^= 0. In order to 
establish the most general formulae, we need two other independent variables, 
t and y. The gradient of expansion variables \ an d ip 1S (auxiliary notation is 



uj = X : ) : 

X, = 1 + fx 2 (6-40) 

Xt - -C + C xX -\(CS + C xx ) X 2 (6.41) 

Xy = -K + K xX -\(KS + K xx )x 2 (6.42) 

(Log^) x = x.- 1 (6.43) 

(LogVOt = -C X - 1 + \c x = -C(Logi>) x + ^C x (6.44) 

(LogVO a = -Kx' 1 + \k x = -A-(LogV)* + \k x (6.45) 

co x - -^ 2 -f (6-46) 
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io t = Clo 2 -C x u + ^(CS + C xx ) = (-Cuj + ^C x ) x (6.47) 
u y = Ku 2 -K x l;+^{KS + K xx ) = (-Klu+^K x ) x (6.48) 



(note that eq. (6.41) generates the eight others) where S,C,K are elementary 



homographic differential invariants linked by the cross-derivative conditions : 

V X X ((<p XX x)t - (<Pt)xxx) = S t + C xxx + 2C X S + CS X = (6.49) 
fx 1 ((<Pxxx)y ~ {<Py)xxx) = S y + K xxx + 2K X S + KS X = (6.50) 
<Px l ((<Py)t ~ ift)v) = C y -K t + C X K - CK X = 0. (6.51) 

Kruskal's choice is implemented by putting S = 0, C = ft, K = f y 



eq. (6.40)-(6.48), thus reducing each rhs to one term and making eq. (6.49)- 
( |6.51 ) useless. 

The function ip — ip>o never appears in the above formulae. Similarly, the 
explicit expressions of Xi i>i C, K as functions of ip — ip® are not needed during 
the computations. We recall them here only for reference : 

X = 1 (6.52) 



_ i 

ip = (ip- tp )tp x 2 (6.53) 



•Pxxx 3 / ip xx \ / tp xx \ 1 / ip 



o r 1 viii " / rii \ / vii \ ^ I fxx i /„ - .\ 

S = {ip; x} = = - - (6.54) 

fx 2 \(p x J \ (p x J x 2 \ ip x J 

<Pxx\ ,f^_ (6 . 55) 



ip — ipu 2ip x J x \(p-(po 2<p 
C = (6.56) 

<Px 

K = -^L. (6.57) 

<Px 

In some applications, it is necessary to choose for x the most general homo- 
graphic transform of fl6.52| ) which vanishes as tp — ipo 

gradx = X +X lX + X 2X 2 , (6.58) 
gradw = -X 2 - Xjw - X w 2 (6.59) 

gradLog^ = XoX^ + ^Xi. (6.60) 

The vectorial coefficients Xj depend on (S, C,K, . . .) and two additional arbi- 
trary functions. The auxiliary expansion variable ip is defined by its logarithmic 
gradient and by the condition that it should vanish as <p — ipo . 
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6.6 The Painleve test 



The synthesis of the different methods to generate necessary conditions for 
the PP produces the following algorithm, called "Painleve test" . 



C onsi der a D E ( p.30|) of order N, already transformed if necessary, see sec- 
tions |6.l| and |6.2| , so as to be polynomial in u and its derivatives, analytic in x. 
The Painleve test is made of the following steps. 



Step 0. Perform a transformation (3.5) in order to reduce the number of terms 
in the equation (details in section |5.3.1 and ref. M.I and, for PDEs, |3j|). 



Ex. : equation u x +u t +u xxt +u x u t = 0, under the translation u = U—x—t, 
becomes U xxt + U x U t + 1 = 0. 

Step 1. Require the satisfaction of the very general necessary conditions ob- 
tained by Painleve (details in section 5.5.2| ). 



Ex. @ : -3u 2 u'u"' + hu 2 u" 2 - uu' 2 u" - u' 4 = 0. Unstable for 5/3 has 
not the required value 1 — 1/n. 

Ex. : (l + u 2 )u xx — 2uu 2 +u 2 = 0. The ODE obtained by the reduction 
(x, t) —> x—ct has an A with two simple poles u — ±i and residues l±?c 2 /2. 
The ODE is unstable, and so is the PDE. 

Step 2. If the degree is greater than one, establish the ODE satisfied by the 



singular solutions (details in section 5.1) 



Step 3. Put the DE under a canonical form of Cauchy; find all the excep- 
tional points where the Cauchy theorem fails; for each such point, define a 
homographic transf orm atio n fl3.5|) allowing the Cauchy theorem to apply 
(details in sections |5.3| and |5.5| ). For each DE (the original one and all 
these homographic transforms), perform step 4. 

Ex. : ( P5) h as the exceptional points u — 1 and u — (poles of A, see 



section 



5.5.2)) . These points are regular for the ODE in (u— 1) 1 and u 1 . 



Ex. : the reduced three-wave interaction dynamical system 

x' = —2y 2 + z + -fx + Sy, y' = 2xy + jy — Sx, z' = —2xz — 2z 

has the exceptional point y = 5/2 El], not so evident on the system itself 
but easily unveiled by considering the equivalent third order ODE for y(t). 

Step 4. Find all the families u ~ u[ ) ' ) x p (uq -* ^ 0) (details in section |5.4| ). 
Discard those families which are also families of the ODE for singular 
solutions established at step 2. Require all components of remaining p's 
to be integer. Discard all families having all components of p positive. 
For each remaining family, perform step 5 and at least one of the steps 6 
and 7. 
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Ex. : (P5) has six families of movable simple pole-like singularities 

(6.61) 

Ex. : the reduced three-wave interaction has the families |Tl| 

(x,y,z) ~ (-(l/2)x _1 ,(«/ 2 )x _1 )«o), zo arbitrary, indices (-1,0,2) 
(ar,tf s «)~(x" 1 »<y/2,-X" 2 ). indices (-1,2,2). 

The first one will pass the test while the second will generate at index 2 
the conditions jS — 0, 7(7 + 1) = 0. 

In case the DE has too many terms, this step is worth being programmed 
on a computer, by fear of missing some families. 

Warning. If one is unsure about some component u of u behaving like a 
positive integer power p of x, it may be safer to switch to the DE for u^ 1 . 

Step 5. From the auxiliary equation of the simplified equation, compute the 
linear operator P(i ) eq . ( [2.36 ) and the indicial equation ( [2.36 ) det P(i) = 



(details in section 5.4). Compute its zeroes (the Fuchs indices). Require 



each index to be integer (details in section 5.11) and to satisfy the rank 



condition (5.113) 



Ex. (|M example 5.B). These are two coupled PDEs with a single family 
whose linear operator P(i) is 

The indices are the zeroes of its determinant (—2,-2,-1,1). For the 
double index i = —2, the rank of P(i) is one, so the system of PDEs is 
unstable. 

2 



Ex. |22|, |15|, ^2| : the equation u xxx — luu xx + llu x = has only one family 
with three indices : p = — 1,Mq ^ = —2, indices (—6,-1,-1). The double 
index — 1 immediately proves the instability. 

Step 6. (NonFuchsian case). If the degree of the indicial polynomial is strictly 
lower than N, and if a particular solution is known in closed form, apply 
the NonFuchsian perturbative method (details in section p. 



Step 7. (Fuchsian case) . Denote p the smallest integer Fuchs index, lower than 
or equal to —1. Define two positive integer upper bounds fcmax and "-max 
representing the cost of the computation to come, see advice below. Solve 
the linear algebraic system (5.118) in the unknown u'™\ (J, n) ^ (0,0), 



for the successive values k = to fc m ax, n = to n m ax with j = k + np: 
whenever j is an index i of multiplicity mult(i), 
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- require the orthogonality condition ( 5.119Q to be satisfied for any value 



of the previously introduced arbitrary coefficients, 

- if n = or (n = 1 and i < 0), assign arbitrary values to mult(i) 

components of defining a basis of Ker P(i), 

- if [n — 1 and i > 0) or n > 2, assign the value to mult(z) components 

of defining a basis of Ker P(i). 



Details in section 5.7 



Advice for choosing k m ax an d ^max : if the order n = fails to describe 
the general solution, take at least nmax = 2; take fcmax so as to test the 
greatest Fuchs index for n — rimax (all details in the remarks at the end 



of section 5.7) 



This ends the test. Step 6 has been put before step 7 because in all our 
examples it allows to conclude sooner. 

Let us again stress that these sets of conditions may not be sufficient : Pain- 
leve gave the counterexample of the second order ODE whose general solution 
is ± sn[ALog(cia; + C2);fc)], with (c\,C2) arbitrary, for which no local test can 
generate the necessary and sufficient stability condition that 2ni\ be a period 
of the elliptic function sn. For advanced features, see section 5.£ and El). 



6.7 The partial Painleve test 



In the search for the tiniest piece of integrability, the physicist, see section 
■7j will perform the above Painleve test to its end, i. e. without stopping even in 
case of failure of some condition, so as to collect a bunch of necessary conditions. 

Turning to sufficiency, these conditions will then be examined separately in 
the hope of finding some global element of integrability, most often a Darboux 
eigenvector. 



For instance, the Lorenz model (1.2) has two families 



x ~ 2i X -\ y - -{2i/a) X ' 2 , z ~ -(2/a) X - 2 , z 2 = -1, (6.63) 
with the same indices (—1,2,4), which generate the no-log conditions [1061 B2) 



Q 2 = (8/3)(6 - 2<r)(b + 3<t - 1) = 

Q 4 = -4i(6 -a- 1)(6 - 6a + 2)x 2 + (8/3) (6 - 1)(6 - 3cr + 1)5 
-Aba(b - 3cr + 5)r + f{b, a) = 0, 

in which x 2 is arbitrary, S is the Schwarzian of the invariant analysis, and / a 
polynomial irrelevant for what we want to emphasize. Performing a logical or 
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operation on these conditions instead of the logical and of the mathematician, 
one obtains the condition on (b, a) 

(b - 2a){b + 3cr - l)(b - a - l)(b - 6cr + 2)(b - l)(b - 3<r + 1) = (B.64) 

What is remarkable is that known analytic results on this model (first in- 
tegrals |7q| , particular solutions [p2| , Darboux eigenvectors [fl8j 0) belong to 
one of these six cases. Conversely, to each of the six factors there corresponds 
such a result, although sometimes only for a finite set of values of (b, a). 
Remarks. 

1. With the restriction S = one would miss two of the six factors. 

2. First integrals P(x,y,x)e xt , with P polynomial and A constant, should 
not be searched for with the assumption P the most general polynomial 
in three variables. Indeed, P must be an entire function of t i. e. have no 
singularities at a finite distance. The generating function of such polyno- 
mials is built from the singularity degrees of (x, y, z) |H] 



(1 — ax)(l — a 2 y)(l — cpz) 

and it provides the basis, ordered by singularity degrees 

(1), (x), (x 2 ,y,z), (x 3 ,xy,xz), (x i ,x 2 y,x 2 z,yz, z 2 ,y 2 ), .(6.66) 

Thus, Pi should be searched for as a linear combination of (1, x, x 2 , y, z). 
All known first integrals are found at the P4 level . 

3. The case b = 1 — 3er is on an equal footing with the case b = 2a which 
admits the first integral (x 2 — 2az)e 2crt , but finding its first integral is still 
an open problem. 
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Chapter 7 

Sufficiency : explicit 
integration methods 



We review the algorithmic methods which may perform the explicit integra- 
tion, with emphasis on ODEs. The PDE case is handled in another part of this 
volume Q. 

We assume that the application of the Painleve test (necessary conditions 
for the PP) has led either to no failure or to a minor failure, corresponding 
respectively to a presumption of integrability in the Painleve sense or of partial 
integrability. If perturbative methods have been used, one has to decide to give 
up at some perturbation order n (remember the counterexample of Painleve). 
The goal is then either to prove the sufficiency (integrability) or to build par- 
ticular solutions (partial integrability). 

If the DE belongs to one of the fully studied classes enumerated in chapter 
W, the question is solved. Indeed, either it is possible, by some homographic 



transformation (3.5), to bring the DE back to a normalized ("classified") DE, 
in which case the integration is finished, or this is impossible, in which case the 
DE has not the PP. 

For a DE which has not been classified, if one excludes the case where the 
DE is an ODE and defines a new function (a quite improbable event which has 
not occured since 1906), the explicit proof of sufficiency amounts to (the cases 
below are not mutually exclusive) 

• either (ODE case) express the general solution as a finite expression of a 
finite number of elementary functions (solutions of linear equations, the 
Weierstrass p function, the six Painleve functions), 

• or (PDE case) find a Lax pair. 

In the partial integrability situation, one tries to obtain degeneracies of these 
results : a particular solution or a pair of linear operators able to generate a 
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subclass of solutions. 

The methods to handle both cases are the same, and they again only rely 
on the singularity structure. Their basic common idea is that the singular part 
of the Laurent expansions (of a local nature) contains all the information for a 
global knowledge of the solution. 

The two existing methods are known as the singular part transformation and 
the truncation method. Before describing them, let us give a few definitions and 
explain how Painleve proved the sufficiency for the six equations (P1)-(P6). 



7.1 Sufficiency for the six Painleve equations 

Painleve introduced the concept of "integration parfaite" and used it to solve 
the question of sufficiency for the six equations discovered by himself and Gam- 
bier. The idea is to perform a finite (in the sense of Poincare : finite expression) 
single valued transformation from (Pn) to another ODE which has no more 
movable singularities although it may still have fixed critical singularities. Such 
an ODE has qualitatively the same singularities than a linear ODE, and Painle- 
ve says that its integration is then "parfaite" (achieved) (BSMF p. 205) : given 
any initial conditions, its solution can be computed with an arbitrary accuracy 
(by e.g. the sequence of coefficients of convergent Taylor series) since one knows 
in advance where the remaining (fixed) singularities arc located. The movable 
singularities of the original ODE are then totally under control. The equations 
with fixed critical points therefore constitute a natural extension to the linear 
equations. 

Painleve defined such transformations (nowadays called "singular part trans- 
formations") for each of the six equations (P1)-(P6). These transformations, via 
logarithmic derivatives, transform (P1)-(P6) into equations for ip without mov- 
able singularities ((PI) Acta p. 14, (P2) Acta p. 15, (P3) Acta p. 16, (P4,P5,P6) 
CRAS 1906, Oeuvres III p. 120) 

(PI) u=-d*Logip (7.1) 

(P2) u = d x Log - d x Log t/> 2 (7.2) 

(P3) u = e- x {d x \,ogi)i-d x Logil)2) (7.3) 

(P4) u = d x Log i> x - d x Log V> 2 (7.4) 

(P5) u = xe- x (2a)- 1 / 2 (d x Log^ 1 -d x Log^ 2 ). (7.5) 

(P6) u = x(x- l)e- x (2a)- 1 / 2 (d x Logi; 1 ~d x Log^ 2 )- (7.6) 

The Lax pairs of (P1)-(P6) can be found in Re f. |p£ 
The two methods developed in next sections ( |7.3| ) and ( |7.4| ) rely on this 
result. 

The logarithmic derivative plays a privileged role, as generator of a movable 
simple pole with a residue generically unity. A prerequisite to the algorithmic 
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derivation of a transformation from u to ip such as (7.1) is the introduction of 
a free gauge function which we denote <p. 

Such a gauge naturally arises if one thinks of an ODE as the canonical 
reduction of a PDE defined by suppressing the dependence upon all independent 
variables but x. This is the function <p used in the description of the movable 
singularities by ( 6.23| ) rather than x — Xq — 0. Useless at the stage of building 
necessary conditions (the Painleve test), this feature is the key to the algorithmic 
explicit integration methods. 



7.2 The singular part(s) 



Definition. The singular part of one of the families of movable singularities 
of a given DE is the finite sum of the Laurent series restricted to the nonpositive 
powers in the method of pole-like expansions 



Ut 



3=0 



j+P 



(7.7) 



Synonyms are : truncation, truncated expansion. 

Given tp, the singular part ut is a one-parameter (<po) family of expressions 
wt(¥>o)j and the two particular values ipa = and ^0=00 are of special interest. 
For the example of KdV (6.26) 

(n s Jvi tpxx] 2 <pt 2 / ip xxx 3 r^i^ ,. s , 

ut(0) = -2 — + (7.8 

ip 2(p x j 6tp x 3 V ip x 2 [ fx J I 

t \ A ^xx] 2 .ft 2 / ip xxx 3 r^i 2 ^ , . 

u T (oo) = -2 - - — + . (7.9) 

2tp x ] 6tp x 3 V ip x 2 [ <p x J J 

Definition. The singular part operator I? of a family is defined by 

Log^3 — * VLogip — ut(0) — ut(oo). (7-10) 

Example 1 (KdV). The operator D is linear and equal to 2d x . This linearity 
is strongly linked with the Darboux transformation [ p4[ . 

Example 2. The single family of (PI) and the two families of (P2) have the 
singular parts 



(PI): u T 



s 



-dl 



(P2) : u T = ± X -\ V = ±d x 



(7.11) 
(7.12) 
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7.3 Method of the singular part transformation 



This is the method used by Painleve and outlined in previous section 7.1 
It consists of transforming the DE for u into a DE for tp by the nonlinear 
transformation 

u = VLogp, (7.13) 

where T> is the singular part operator associated to one of the families of the 
equation for u. 

If the transformed equation for tp. can be integrated, so is the original equa- 
tion. 

Example 1 (linearization) . The unique first order first degree ODE with the 
PP, namely the Riccati equation ( 1 



LI), h as a V oper ator e qual to — a 2 d x , com- 



putable from the basic formulae (3.40), ( |6.52 ) and ( 7.1 0| ) . The transformation 



hog tp from u to tp leads to the second order linear equation (1.2) 
for tp. It is then sufficient to know two particular solutions tpi and tp 2 (which 
are functions) of this linear equation to have a global knowledge of the general 
solution of the Riccati equation by the formula 

u = -a 2 1 d x Log(ci<,5i + c 2 tp 2 ). ( 7 -!4) 

Similarly, the transformation p = —d 2 Log a associates to the Weierstrass 
elliptic function p a function a which is an entire function, solution of a nonlinear 
ODE. 



Example 2 (simplified equation of one of the 50 stable ODEs (4.8)). The 
ODE 

E = u" + uu'-u 3 = (7.15) 

possesses two families of movable simple poles uq = 1 and uo — —2, with the 
one-parameter particular solutions uq/{x — xq). The first family operator is 
T> = d x and it transforms it into 

u = 9 x Log^ E = tp(^?pj =0, (7.16) 

which integrates as tp = ap(x — xq, 0, 173) with (a, xq, 173) arbitrary and provides 
the general solution. The two families of movable simple poles for u correspond 
to the movable simple zeroes of p (residue uq = 1) and to the movable double 
poles of p (residue uo = —2). 

Example 3 (indirect linearization). The Ermakov-Pinney equation ( |p.4| ). 
after the transformation u~ 2 = v removing its algebraic singularity 

E=-^vv xx + ^v 2 x ~a 2 v 2 +{3 2 v 4 = 0, (7.17) 
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has two families v ~ ±(2/3) 1 x 1 , and the transformed ODE under 
v = (2p)- 1 d x Logi P M 



Wxxx 3 ( \p 



fx 2 V V: 



xx 1 o„2 



= -2a 2 (7.18) 



is a Schwarz ODE fl6.6| ). This integrates the Ermakov-Pinney equation via a 
finite two-valued expression. 

Well suited to DEs possessing only one family (Riccati, Weierstrass, (PI), 
KdV), this transfor matio n must be adapted, following the Painleve formulae for 



(P2)-(P6) in section |7.l| , to suit DEs with more than one family (Jacobi elliptic 



equation, (P2) to (P6)). This is done in the course on PDEs 



7.4 Method of truncation (Darboux transforma- 
tion) 

Perfectly adapted to PDEs [Q, this method is rather poor for ODEs, for 
an intrinsic reason which is the absence in this case of a Backlund transforma- 
tion (link between two different solutions of the same DE introducing at each 
iteration at least one more arbitrary parameter in the solution) . It nevertheless 
succeeds, at least partially, in many situations. 



The idea [ 1 1 1 1 is to consider the singular part (7.7) of one family (or the 
sum u = J2f D Log ipf of the singular parts of several families) as a parametric 
representation of a solution in terms of one function ip linked to \ by X -1 = 
d x Log^? (or several functions tpf, one per family /). Every function ipf, which 
defines a singular manifold ipf = 0, is required to be an entire function, and for 
instance to satisfy the same linear system of two PDEs L\ — 0, L2 — with 
some adjustable coefficients. 

The method consists of identifying to zero the lhs E(ut) considered as a 
polynomial of ipf and its independent derivatives modulo the constraint that 
each ipf satisfies the linear system. This generates an overdetermined set of 
determining equations whose unknowns are the coefficients Uj of ( [7.7[ ) and the 
coefficients of the linear system. The remarkable fact is that the determining 
equations are easy to solve. 

The result is some class of exact solutions, and this class is easily interpreted. 
If the commutator [L\,L2\ is identically zero (which is always the case if the 
linear system has constant coefficients), the solutions are particular ones (PDE 
case) or any kind (particular or general) (ODE case). If this commutator is 
zero only when some coefficient of {L\,L2) satisfies some PDE, quite probably 
(Li, L2) define a Lax pair. 

Again, the ODE case to which we restrict is much less rich than the PDE 
case [p4| to which we refer the reader. 
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7.4.1 One-family truncation 



This is the celebrated WTC truncation procedure [ 1 1 1 1 . Applicable to any 
DE with any number of families, it consists of selecting one of the families ip = 0, 
in which tp obeys the linear system of the invariant analysis 

^+2^ = 0, (7.19) 

A + C^ x - ^ = 0. (7.20) 

The functions S et C are adjustable fonctions, only constrained by the cross- 
derivative condition ( 6.49| ). Consider for instance the Ermakov-Pinney ODE 

E = -I Wxx + ^-aV+/?V = 0. (7.21) 

The infinite Laurent series is v = (2/3) % +«i + 0{x) with v\ arbitrary and 
(3 one of the two square roots of (3 2 . Thanks to the gauge (p, the coefficient v\ 
is not a constant but a function. 

The method consists of assuming that a solution v can be represented by 
the truncation 

v = v T = T^X" 1 +vi =VLogiP + v 1 (7.22) 
implying for the lhs E of the DE the similar truncated expansion 

4 

E=J2E jX j -*- (7-23) 

i=o 

This generates, in this example, five equations Ej = in the unknowns (v\, S). 
Among them, Eq is zero since the coefficient vq of the series for v is already the 
good one. E\ is zero since 1 is a Fuchs index whose orthogonality condition is 
satisfied. Denoting vi = —Vi/(20), there remain the three equations 

16(3 2 E 2 ee -4a 2 + S + W 2 + 6Vi, x = 0, (7.24) 
16P 2 E 3 ee 8a 2 Vi + 2SVi - 4Vf + S x + 2V 1<XX = 0, (7.25) 

16/3 2 E 4 ee ^S 2 - 4a 2 V? + - V X S X + 3SV 1>X + 3V? jX - 2V 1 V 1 , XX t=?H6) 

The algebraic elimination (i. e. without differentiation) of V\. x and Vx lXX 
among these three equations yields (S — s) 2 = 0, with s = —2a 2 , then V\ is 
found to satisfy the Riccati equation 

-2V hx - 2V? = s. (7.27) 
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Hence the particular solution 

v = 



-Vi), 



(7.28) 



in which each variable x _1 and V\ satisfies the same Riccati equation and de- 
pend on one arbitrary parameter. This is the general solution, which can be 
written as v — (2f3)~ 1 (d x Log-0i — c^LogV^) in agreement with the structure 
of singularities, cf. (\l.2\). 

Remark. The class of particular solutions generically found by this method 
is the class of polynomials in tanh, which correspond to a constant value for S. 

Another example is the (P2) equation, for which the one-family truncation 
u = x^ 1 + ui provides the one-parameter particular solution u\ = 0, S — x on 
the condition a = 1/2, i. e. an algebraic transform of the Airy equation. 



7.4.2 Two-family truncation 



When a DE admits two families with opposite principal parts, such as (7.17), 
it is natural to seek particular solutions described by two singular manifolds J34| 

v = -!-r[e?xLogV>i - SxLog^a +«o], (7.29) 
Zp 

in which (ipi,ip2) is a basis of the two-dimensional space of solutions of some 
ODE whose general solution is entire, e.g. the second order linear equation with 
constant coefficients 

k 2 

- -jV = (7.30) 

* 2 = Cie** + C 2 e~i' = C cosh^(x-a;o), (7.31) 
ipi(x) = 1 $2(x + a), i>i{x) = ^2{x — a), a arbitrary, (7.32) 



Substituting (7.29) into (7.17) and eliminating any derivative of (V ; i;V ; 2) of 
order higher than or equal to two in x results into a polynomial in the two 
variables ipi^x/ipii^.x/^- Before identifying it to the null polynomial, one 
must take account of the first integral /io, the ratio of two constant Wronskians 

1pl,xi>2,x k 2 k ( tp hx 1p2,x\ , , , ,„ QQ x 

M07: I —r f— , Co = cotanh/ca, (7.33) 



ipi ip 2 4 2 \ Vi i>\ 

which splits the polynomial of two variables into the sum of two polynomials in 
one variable : 



16/3' E = (k 2 - 4a z + 6v^ + 6kfj, v )( 



2 / , \ 2 



1pl,x\ . ( fp2,x 



i>l J \ i/>: 
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+(kno(k z - Aa z + &Vq) + 2(3fc>o - k - 4a 2 + 4v£)v ) 



Vi 



¥>2 



+2a 2 k 2 - — - 3fe 3 MoW - 4a 2 u 2 - 3fc 2 v 2 + v$. (7.34) 

This defines three different algebraic equations in the unknowns (k,vo, fio); 
their two solutions 



k 2 
k 2 



Aa , vq = 0, fiQ arbitrary, 
4a 2 , vq = 2a, fc/io = —2a 



(7.35) 
(7.36) 



are just two different representations |3J] of a solution of ( 7.17 ) depending on 
two arbitrary constants (fj,o, xq) ■ with this simple assumption, we have obtained 
the general solution 



1 

2(3 



Ipl 4>2 



sinh kc 



(3 cosh fc (x — xq) + cosh ka 



(7.37) 



In particular, with fiQ = one thus obtains immediately the class of solutions 
polynomial in the two variables tanh and sech |33|, thus augmenting the class 
indicated at the end of previous section. Evidently, if the DE has only one 
family, no dependence on sech can be found. 
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Chapter 8 

Conclusion 



The solution of an ODE cannot escape the structure of singularities of the 
ODE. Such a structure can be studied on the equation itself, without any a 
priori knowledge of the solution, providing a deep insight on the possibility or 
not to perform the explicit integration. 

Two levels of integrability have been defined : the Painleve property (the 
most elementary level) and the integrability in the sense of Poincare (the prac- 
tical level). 

A first series of methods (globally called "the Painleve test" ) provide neces- 
sary conditions for a differential equation to have the Painleve property, without 
any guarantee on the sufficiency. In case of a negative answer from these first 
methods, there exist other methods (Legons 8, 9, 10, 13, 19), not developed 
here, to provide necessary conditions for the general solution to have only a 
finite amount of movable branching, which implies the integrability in the sense 
of Poincare, a weaker property than the PP. 

In case of a positive answer, the DE may have the PP, i.e. a general solu- 
tion free from movable critical singularities. Then, a second series of methods 
are available to perhaps constructively prove the PP by explicitly building the 
general solution or some equivalent information (Lax pair). In case of failure of 
these second methods, the only remaining tool is human ability. 

There exists another approach to DEs which is not based on the study of 
singularities, this is the method of infinitesimal symmetries |87| , |88|| . It provides 
reductions of PDEs to "smaller" PDEs or to ODEs, and it may provide first in- 
tegrals of ODEs. However, the PDEs or ODEs left over after its completion still 
require to be integrated, and the only methods to do so are those based on sin- 
gularities. For instance, with the ODE (PI), the method of symmetries cannot 
provide any information (existence or not of a first integral, single valuedness 
or multivaluedness). 
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